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Trial  runs  of  a  1:24  scale  model  crane  ship  were  conducted  in  the  David  Taylor  Model 
Basin.  The  model’s  response  to  regular  waves  under  various  ship  configurations,  crane 
configurations,  sea  states  and  ship  headings  relative  to  the  incoming  waves  were  recorded.  The 
Response  Amplitude  Operator  (RAO)  Program  analyzes  the  frequency  responses  to  controlled, 
regular  waves  and  generates  full-scale  RAOs  as  a  prediction  of  the  actual  ships  response. 
Accurate  generation  of  these  full-scale  RAOs  enables  future  prediction,  using  the  principle  of 
linear  superposition,  of  ship  motions  in  an  irregular  sea  to  be  compared  to  actual,  full-scale  trial 
runs  being  conducted  off  the  coast  of  California  near  Camp  Pendleton  in  September  2000. 


v 


TABLE  OF  CONTENTS 


I.  INTRODUCTION 


1 


II.  PROBLEM  FORMULATION . 5 

A.  THEORY . 5 

1.  Random  Process . . 5 

2.  Energy  Spectra  and  Power  Spectral  Densities. . 5 

3.  Response  Amplitude  Operators . . 6 

B.  COMPARISON  TOOLS - 7 

1.  Spreadsheet  Results.............................................. . 7 

2.  MATLAB  Methods . . 8 

a.  Welch . 8 

b.  Periodogram . 10 

c.  Multi-taper . 11 

III.  PROGRAM  DEVELOPMENT . . 13 


A.  GENERAL  DESCRIPTION . . 13 

1.  Initial  Set-up ...... - .... - .......... - ........ — - -  13 

2.  Establish  Proper  String  Configuration  for  the  Data  Run .......  15 

3.  Determine  Number  of  Data  Channels  to  be  Processed ............  15 

4.  Establish  File-path  and  Load  the  Raw,  Time-based  Data.......  16 

5.  Compute  both  Model-Scale  and  Full-Scale  Omega  Ranges ....  17 

6.  Full-Scale  Scalers,  PSD  Units  and  RAO  Units  Library . 17 

7.  Looping/Processing  of  Data  Channels.......................................  18 

8.  Base-lining  the  Raw,  Time-based  Data . 20 

9.  FFT  Computation - - ... . 21 

10.  Model-Scale  PSD  Computation - ...... - ..... - ............  21 

11.  Moving-window  Averaging  of  Model-Scale  PSDs ....................  22 

12.  Computation  of  Model-Scale  RAOs  — . . 26 

13.  Computation  of  Full-Scale  RAOs ..............................................  27 

14.  Plot  Model-Scale  PSDs  and  Save  Figure  in  jpg  Format .... —  27 

15.  Plot  Model-Scale  RAOs  and  Save  Figure  in  jpg  Format .........  28 

16.  Plot  Full-Scale  RAOs  and  Save  Figures  in  jpg  Format ...........  29 

17.  Save  Omega  Matrices  as  Last  Column  of  Composite 

PSD  and  RAO  Matrices ..............................................................  31 

18.  Save  Composite  PSD  and  RAO  Matrices.................................. 31 


Vll 


B.  VALIDATION _ 32 

1.  Load  Applicable  Matrices............ — ...........................................  32 

2.  Base-line  the  Raw,  Time-Based  Data.........................................  33 

3.  Compute  Wave  Ht  PSDs  using  MATLAB  Methods .......... 33 

4.  Compute  TACS  Roll  PSDs  using  MATLAB  Methods.............  34 

5.  Convert  MATLAB  Method  Frequencies  to  Omegas ... — ........  34 

6.  Compute  Areas  under  the  PSD  Curves  for  Comparison _ 35 

7.  Graphical  Comparison  of  Model-Scale  PSDs ...........................  37 

8.  Compute  Model-Scale  RAOs........................... . 41 

9.  Graphical  Comparison  of  Model-Scale  RAOs... — ...................  42 

10.  Scale  to  Full-Scale  RAOs  and  Omegas ......................................  43 

11.  Graphical  Comparison  of  Full-Scale  RAOs....................... _ ...  44 

IV.  CONCLUSIONS  &  RECOMMENDATIONS _ 47 

A.  CONCLUSIONS _ 47 

B.  RECOMMENDATIONS _ 48 

APPENDIX  A.  TEST  CONFIGURATIONS  I _ 49 

APPENDIX  B.  TEST  CONFIGURATIONS  II _ 61 

APPENDIX  C.  TEST  CONFIGURATIONS  II _ 63 

APPENDIX  D.  MASTER  CHANNEL  LIST:  CONFIGURATION  I _ 65 

APPENDIX  E.  MASTER  CHANNEL  LIST:  CONFIGURATION  II _ 67 

APPENDIX  F.  MASTER  CHANNEL  LIST:  CONFIGURATION  HI _ 69 

APPENDIX  G.  RAO  PROGRAM _ 71 

APPENDIX  H.  PTAVG  COMPARISON  PROGRAM _ 89 

APPENDIX  I.  PSD  AND  RAO  COMPARISONS  PROGRAM. _ i, _ 91 

LIST  OF  REFERENCES _ 97 


vm 


BIBLIOGRAPHY... . 99 

INITIAL  DISTRIBUTION  LIST . 101 


IX 


I.  INTRODUCTION 


From  July  15th  to  August  17th  of  1997,  trial  runs  using  1:24  scale  models  of  a  T- 
ACS  Auxiliary  Crane  Ship,  a  DDG-963  Class  Ship,  a  Commercial  Container  Ship  and 
Lighter  Barges  were  conducted  in  the  seakeeping  basin  of  the  David  Taylor  Model  Basin 
at  Carderock  Division  Naval  Surface  Warfare  Center  (CDNSWC),  Bethesda,  Maryland. 
T-ACS  Auxiliary  Crane  ships  are  ran  by  the  Military  Sealift  Command  Ready  Reserve 
Force.  These  self-sustaining,  rapidly-deployable  ships  support  military  sea  transportation 
needs  and  are  extremely  useful  in  ports  that  have  limited,  damaged  or  undeveloped  port 
facilities.  The  functions  of  T-ACS  ships  are  to  lift  and  transfer  various  loads  from  either 
themselves  or  adjacent  vessels  and  piers. 

A  total  of  251  different  model-scale  trial  runs  were  performed,  each  running 
approximately  8  minutes,  under  a  matrix  of  variable  conditions: 


Variable  Condition: _ Description: _ 

Ship  Configurations  Config  I  T-ACS  in  Center 

Container  Ship  to  Port 
Lighter  Barges  to  Starboard 
Config  II  T-ACS  in  Center 

DDG-963  Class  Ship  to  Port 
Lighter  Barges  to  Starboard 
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Config  m  T-ACS  in  Center 


Ship  Heading 

Sea  State 

Boom  Slew  Angle 
Rider  Block  Location 
Rider  Block  Inhaul  Angle 
Boom  Luff  Angle 


Lighter  Barges  to  Port 

45  degree  increments  from  0  to  360  degrees  relative 
to  incoming  waves 

Model-scale  sea  states  3,  3  +  swell,  4  and  4  +  swell 
Angles  of  0, 45,  90,  270  and  315  degrees 
Various  positions  from  0  to  45.3  degrees 
Various  angles  from  0  to  16  degrees 
Angles  of  25,  29.6,  54.5  and  60  degrees 


A  composition  of  wave  height,  body  motion,  velocity  and  acceleration  data  for  the 
vessels  was  recorded  for  each  run  from  an  array  of  sensors  with  their  locations 
dependant  upon  the  ship  configurations  listed  above.  Results  from  the  213 
Configuration  I  trial  runs,  the  most  probable  ship  configuration,  were  provided  for 
frequency  response  analysis  of  the  raw  data  to  estimate  the  Full-Scale  Response 
Amplitude  Operators  (RAOs)  of  the  T-ACS  auxiliary  crane  ship.  An  RAO  is  basically  a 
measure  or  ratio  of  a  vessels  response  to  a  regular  wave  of  unit  amplitude  and  thus 
defined  accordingly: 


amplitude  _of  _  response 
amplitude  _of  _  the  _  wave 


(Zubaly,1996,pp.  322) 


2 


Tupper  expresses  the  method  of  studying  responses  in  a  seaway  in  the  following 


manner: 


"...  the  apparently  random  surface  of  the  sea  can  be  represented  by  the 
summation  of  a  large  number  of  regular  sinusoidal  waves,  each  with  its 
own  length,  height,  direction  and  phase.  ” 

(Tupper,  1996,pp.  104) 

"...  the  response  of  the  ship  in  such  a  sea  could  be  taken  as  the  summation 
of  its  responses  to  all  the  individual  wave  components.  Hence  the  basic 
building  block  for  the  general  study  of  motions  in  a  seaway  is  the  response 
to  a  regular  simusoidal  wave.  ” 

(Tupper,  1996,pp.  104) 


These  concepts  strongly  identify  the  need  for  multiple  trial  runs  using  a  matrix  of  variable 
conditions.  The  study  of  responses  such  as  severe  rolling  in  a  beam  sea  or  excessive 
pitching  and  heaving  in  a  head  sea  provides  extremely  valuable  insight  into  a  vessels 
limitations.  Operating  or  load-handling  contraints  can  then  be  established  based  upon 
existing  as  well  as  predicted  conditions  to  prevent  material  damage,  downtime  and/or, 
even  worse,  human  injury. 
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II.  PROBLEM  FORMULATION 


A.  THEORY 


1.  Random  Process 

Determining  RAOs  begins  with  periodic  sampling  of  random  processes  such  as  a 
changing  waveheight  or  a  vessels  changing  responses  to  the  changing  waveheight  over  an 
interval  of  time.  In  this  case,  each  data  channel  was  sampled  at  a  rate  of  32.2  Hz  over  an 
interval  of  approximately  8  minutes  per  data  run.  This  sampling  produces  the  raw  data 
from  which  the  entire  RAO  determination  process  is  based  upon.  But,  before  the  raw 
data  can  be  used,  it  must  be  base-lined  to  remove  any  offset  or  bias  errors  inherent  in  the 
sampling  system  to  prevent  them  from  contaminating  the  process.  This  is  done  by 
determining  the  mean  of  the  data  points  for  a  given  channel  and  then  subtracting  this 
mean  from  each  data  point.  When  this  step  is  completed  for  each  data  channel,  the 
process  is  ready  to  continue. 


2.  Energy  Spectra  and  Power  Spectral  Densities  (PSDs) 

The  next  phase  in  determining  RAOs  is  to  generate  a  representative  energy 
spectrum  curve  for  each  data  channel  comprised  of  Power  Spectral  Densities  (PSDs)  vs. 
wave  frequencies  (omegas)  in  which  the  area  under  the  curve  represents  the  energy 
associated  with  the  changing  data.  (Papoulias, 2000, pp.  109)  The  PSDs  for  this  curve  are 
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derived  by  expansion  of  the  data  points  per  channel  into  a  Fourier  series  which 
approximates  the  shape  of  the  energy  spectrum  at  the  various  omegas. 
(Papoulias,2000,pp.  121)  The  mean  square  of  the  base-lined,  raw  data  record  within  a 
narrow  band  of  omega  centerd  at  omega  is  represented  as  follows: 


r  (co)  =  S{(0)h.C0  where  r  =  data  _  mean  _over_Q)_  band 

S  (of)  =  energy  _  spectra 
A  CO  =  omega  _band 
(Papoulias,2000,pp.  109) 


Thus,  as  an  integrity  check  along  the  way  which  is  utilized  later  in  the  PSD  and  RAO 
comparison  porgram,  the  mean  square  of  the  whole,  base-lined,  raw  data  record  should 
equal  the  area  under  the  entire  spectral  curve: 


r  (a>)  =  j S(o))da>  (Papoulias,2000,pp.  109) 

o 


3.  Response  Amplitude  Operators  (RAOs) 

Finally,  the  RAOs  are  determined  by  the  ratio  of  the  response  PSDs  over  the 
reference  waveheight  PSDs  as  follows: 
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where 


RAOr  =  response  _  RAO 


RAOr  |2 


PSD, 

PSDW 


PSDr  =  response  _  PSD 
PSDW  =  waveheight  _  PSD 
(Papoulias,2000,pp.  124) 


B.  COMPARISON  TOOLS 


1.  Spreadsheet  Results 

Spreadsheet  Results  are  selected  PSDs,  RAOs  and  Omegas  for  specific 
Configuration  I  trial  runs  that  were  initially  processed  by  Mr.  Dan  Hayden  of  CDNSWC 
using  a  spreadsheet  program.  Only  seven  data  channels  were  processed  per  run.  These 
channels  were  the  four  pitch  and  roll  channels  (channels  3-6)  as  well  as  the  three  relative 
TACS/Lighter  positions  (channels  28-30).  As  for  the  specific  data  runs,  only  19,  27,  57, 
67,  390  and  392  were  processed. 

The  Spreadsheet  Results  provided  a  reference  to  which  the  RAO  Program  results 
would  be  compared.  If  the  RAO  Program  results  were  consistant  with  each  of  the  seven 
processed  channels,  it  was  assumed  that  the  remaining  channels  would  be  correct  as  well. 
In  order  to  ensure  that  the  referenced  Spreadsheet  Results  were  correct,  three  additional 
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MATLAB  methods  of  computing  PSDs  were  used  to  validate  both  the  Spreadsheet 
Results  and  the  RAO  Program. 


2.  MATLAB  Methods 

In  addition  to  the  FFT  method  in  the  RAO  Program,  three  additional  methods  of 
computing  PSDs  were  utilized  as  comparison  tools  to  ensure  validity  of  both  the 
Spreadsheet  Results  and  the  RAO  Program:  Welch,  Periodogram  and  Multi-taper 
Methods.  The  PSDs  are  computed  in  units  of  power  per  radians  while  the  frequencies  in 
Hz  are  later  converted  to  radians  per  second.  Although  there  are  many  similarities,  each 
method  performs  the  task  quite  differently. 

a.  Welch:  [PSDs,Freqs]  =  p welch ( x,nwindow, no verlap, n fftjs ) 

The  selectable  options  when  utilizing  the  Welch  method  are  as  described 

as  follows: 


Option: _ Description _ 

x  Base-lined,  Raw,  Time-based  Data  Vector 

nwindow  Hamming  Window  Length  for  Modified  Periodograms 
(Default  of  []  divides  data  vector  into  8  equal-length 
windows  with  residual  data  points  discarded) 

noverlap  Number  of  Overlapping,  Windowed  Data  Points  (Default 

of  []  uses  50%  overlap  from  window  to  window) 
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nfft 


Length  of  the  Fourier  Transform  for  each  Window 


fs  Data  Sampling  Frequency  (32.2  Hz) 

(SPT, welch) 


The  Welch  method  produces  PSDs  by  averaging  periodograms  of  overlapping, 
Hamming-windowed  sections  of  the  data  vector.  Hamming  windowing  reduces  possible 
sidelobes  in  the  spectral  estimate  in  order  to  reveal  the  presence  of  weaker  components  of 
the  signal  spectrum  that  may  otherwise  get  hidden.  (Marple,1987,pp.  132)  Specified- 
length,  discrete  Fourier  transforms  (DFTs)  of  the  overlapping  windows  are  computed  as 
follows: 


X  (k  + 1)  =  where  WN  = 

n= 0 

N  =  length(x) 

X  =  DFT 
(SPT, welch) 


Indexing  of  (n+1)  and  (k+1)  is  used  since  MATLAB  vectors  run  from  1  to  N 
instead  of  0  to  N-l.  (SPT,fft)  A  modified  periodogram  for  each  windowed  segment  is 
then  computed: 
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(\  N  « 


Modified? eriodogram  =  S(e1C0)  = 


n 


1=1 


-lZH2 

{nP 


(SPT,periodogram) 


The  individual  periodograms  are  averaged  to  produce  just  one  representative 
periodogram.  (S PT,periodogram )  Finally,  the  PSDs  are  computed: 


PSDs  = 


(\_ 

Js 


W) 

J 


(SPT, periodogram) 


b.  Periodogram:  [PSDs,Freqs]  =  periodogram  (x, window, n fftfs ) 

The  selectable  options  when  utilizing  the  Periodogram  method  are  as 

follows: 


Option: _ Description _ 

x  Base-lined,  Raw,  Time-based  Data  Vector 

window  Window  Coefficients  for  Modified  Periodogram  of  the 

Input  Matrix  (Used  for  “m”x”n”  matricies.  Default  of  [] 
implies  single  column  data  vector) 
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nfft 


Length  of  the  Fourier  Transform 


fs  Data  Sampling  Frequency  (32.2  Hz) 

(SPT,periodogram) 


The  Periodogram  method  produces  PSDs  in  the  same  general  manner  as  the 
Welch  method  previously  discussed.  A  modified  periodogram  is  generated  which  leads 
directly  to  computation  of  the  PSDs.  The  differences  between  the  two  are  that  no 
windowing/segmenting  of  the  data  vector  occurs  and  only  one  periodogram  is  generated. 
(SPT,periodogram)  Thus,  the  Welch  method  is  essentially  a  refinement  of  the 
Periodogram  method. 


c.  Multi-taper:  [PSDs,Freqs]  =  pmtm{x,nw,nfftjs) 

The  selectable  options  when  utilizing  the  the  third.  Multi-taper,  method 

are  as  follows: 


Option: _ Description _ 

x  Base-lined,  Raw,  Time-based  Data  Vector 

nw  Determines  Number  of  Discrete  Prolate  Spheroidal 

Sequences  (n  =  2nw-\)  used  as  Data  Tapers  for  of 
Estimation  of  PSDs  (Default  of  []  is  nw  =  4 ) 
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nfft  Length  of  the  Fourier  Transform 

fs  Data  Sampling  Frequency  (32.2  Hz) 

(SPT,pmtm) 


The  Multi-taper  method  is  by  far  the  most  complex  of  the  three  comparison 
methods  for  production  of  PSDs.  This  method  combines  linear  and  nonlinear  modified 
periodograms  to  estimate  the  PSDs  by  computing  each  periodogram  using  a  sequence  of 
orthogonal  tapers  or  windows  in  the  frequency  domain  as  specified  from  the  discrete 
prolate  spheriodal  sequences.  (SPT,pmtm) 
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III.  PROGRAM  DEVELOPMENT 


A.  GENERAL  DESCRIPTION 

The  RAO  program  contains  a  sequence  of  events  utilizing  standard  MATLAB 
functions.  The  validation  methods  to  be  discussed  later  are  dependant  upon  the 
MATLAB  Signals  Processing  Toolbox  from  which  the  commands  are  comprised  of 
standard  functions  within  the  source  code.  Each  significant  event  of  the  RAO  program  is 
addressed  below: 


1.  Initial  Set-up 

Before  processing  a  data  run,  both  the  raw  data  disc  and  the  output  storage  disc 
must  be  placed  in  their  respective  designated  drives.  Designation  of  these  drives  must 
conform  to  the  system  being  used  for  processing.  In  order  to  start  the  RAO  program,  the 
operator  enters  the  desired  run  number  to  be  processed  as  a  MATLAB  command  in  one 
of  the  following  two  formats: 


Command  Format: _ Applicable  DTMB  Runs: 

rao(##)  19  -  99 

rao  (###)  101-587 
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Once  this  has  been  initiated,  the  program  designates  which  drive  the  raw,  time-based  data 
is  to  be  read  from  as  well  as  which  drive  the  generated  jpg  format  figures  and  composite 
matrices  for  the  specific  run  are  to  be  saved  to.  A  floating-point  format  with  5  digits  is 
placed  into  effect  and  the  following  constants  are  established: 


Constant: _ 

ptavg  =  50 
window  =  500 
lambda  =  24.175 
freq  =  32.2 


Definition: _ 

Number  of  Sequential  PSDs  used  for  Smoothing 
Number  of  PSDs  within  a  Rational  Freq  Range 
Actual  Model  Scale  (1:24.175) 

Raw  Data  Sampling  Frequency  in  Hz 


Selection  of  the  “ptavg”  constant  is  the  unavoidable  compromise  between  display 
resolution  and  accuracy  which  will  be  addressed  later.  As  for  the  “window”  constant,  it 
envelops  the  first  500  smoothed  PSDs  which  encompass  a  reasonable  range  of  omegas  to 
be  addressed  under  full-scale  conditions.  The  remaining  constant,  “freq”,  is  the 
frequency  at  which  the  raw,  time-based  data  was  sampled. 


%  Initial  Set-up - 

function  rao(run) 

sread  =  char ( 'H:  \  ' ) ; 
swri te  =  char( 'D: \ ' 

format  short  e 
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ptavg  =  50; 
window  =  500; 
lambda  =  24.175; 
freg  =  32.2; 

2.  Establish  Proper  String  Configuration  for  the  Data  Run 

Based  upon  the  run  number  entered  by  the  operator,  proper  string  configurations 
for  the  run  as  well  as  the  folder  to  which  results  are  to  be  saved  are  established.  The 
program  then  differentiates  between  a  two  digit  and  a  three  digit  run  number  so  that  the 
run  string  evolves  as  “R###”  and  the  folder  as  “###”.  Broad  applicability  of  designators 
such  as  these  are  utilized  throughout  the  RAO  program. 


%  Establish  Proper  String  Configuration  for  Data  Run — 
if  run  <=  99 

=  char('RO'); 

=  char  ('  O'); 


si 

s2 

else 

si 

end 


=  char  (  'R ' )  ; 


srun 

filename 


num2str (run 
=  strcat (si, srun) ; 


if  run  <=  99 
folder  = 
else 

folder  = 

end 


strcat (s2,  srun) ; 
srun; 


3.  Determine  Number  of  Data  Channels  to  be  Processed 

The  run  number  entered  enables  identification  of  the  number  of  data  channels  to 
be  processed  due  to  the  particular  data  channel  configuration  for  the  run  as  follows: 


15 


DTMB  runs: 


Configuration: 


Data  Channels: 


19-503 

I 

46 

504  -  533 

II 

38 

534-587 

in 

24 

Each  data  channel  is  individually  labeled  per  configuration  as  listed  in  Appendixes  A,  B 
and  C. 


%  Determine  Number  of  Data  Channels  to  be  Processed - 

if  run  <=  503 

channels  =  46; 
elseif  run  <=  533 
channels  =  38; 
else 

channels  =  24; 

end 


4.  Establish  File-path  and  Load  the  Raw,  Time-based  Data 

Given  the  run  number  and  drive  designation  from  which  the  data  is  to  be  read, 

“filepath”  is  established  as  a  string  in  order  to  execute  the  MATLAB  command  “load”  to 

load  the  raw,  time-based  data  matrix  for  processing.  Once  the  matrix  is  loaded,  its 

dimensions  are  identified  utilizing  the  MATLAB  command  “size.”  The  number  of 

columns  representing  the  individual  data  channels  is  designated  as  the  variable 

“channels”  and  the  number  of  rows  representing  the  number  of  data  points  collected  is 

designated  as  the  variable  “n”,  both  of  which  will  be  used  throughout  the  program. 
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%  Establish  Filepath  and  Load  the  Raw,  Time-based  Data 

filepath  =  [sread,  filename]  ; 

load] filepath) 

szdata  =  size (data); 

channels  =  szdata  (2) ; 

n  =  szdata  (1); 


5.  Compute  both  Model-Scale  and  Full-Scale  Omega  Ranges 

Before  going  any  further,  the  Model-Scale  omega  range  is  computed  based  upon 
the  sampling  frequency,  total  number  of  data  points  and  the  “window”  length  designated 
earlier.  Dividing  the  Model-Scale  range  by  the  square  root  of  the  model  scale  “lambda” 
produces  the  Full-Scale  omega  range.  This  is  supported  by  the  basic  assumption  for 
presentation  of  motion  data  that  “Natural  periods  of  motion  vary  as  the  square  root  of  the 
linear  dimension.”  (Tupper,1996,pp.l06)  In  addition,  Tupper  identifies  that  “In  watching 
model  experiments  the  motion  always  seems  rather  ‘rapid’  because  of  the  way  period 
changes  .  Thus,  a  1/25  model  will  pitch  and  heave  in  a  period  only  a  fifth  of  the  full- 
scale  ship. 


%  Compute  both  Model-Scale  and  Full-Scale  Omega  Ranges 
momega  =  2 *pi *freq/n * ( 0 : window- 1) ; 

f omega  =  momega/ sgrt (lambda) ; 


6.  Full-Scale  Scalers,  PSD  Units  and  RAO  Units  Library 

There  are  four  individual  sets  of  Full-Scale  scalers,  PSD  units  and  RAO  units  in 
the  library  associated  with  the  type  of  data  collected  by  the  a  particular  channel.  The 
combinations  are  as  follows: 
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Data  Tvpe: 

Full-Scale  Scaler: 

PSD  Units: 

RAO  Units: 

Wave  Height 

1 

in 

dimensionless 

Distance 

12/lambda 

deg 

deg/in 

Velocity 

12/lambdaA(3/2) 

deg/s 

deg/in  s 

Acceleration 

12/lambdaA(2) 

deg/sA(2) 

deg/in  sA(2) 

%  Full-Scale  Scalers,  PSD  Units  and  RAO  Units  Library- 

fsl  =  1; 

pul  =  '  (in)  '  ; 

rul  =  '(dimensionless)'; 

fs2  =  12 /lambda; 
pu2  =  '  ( deg)  '  ; 

ru2  =  '(deg/in)'; 

fs3  =  12 /lambda*  (3/2)  ; 
pu3  =  '(deg/s)'; 

ru3  =  '  (deg /in  s)  '  ; 

fs4  =  12 /lambda*  2 ; 
pu4  =  ' (deg/s*2)  ' ; 

ru4  =  ' (deg/ in  s*2)  ' ; 


7.  Looping/Processing  of  Data  Channels 


It’s  now  time  to  begin  looping  to  process  the  “channels”  sequentially.  From  the 
first  data  channel  to  the  last,  each  is  assigned  a  Full-Scale  Scaler,  “fscaler”;  a  PSD  Units, 
“psdunits”;  an  RAO  Units,  “raounits”;  and  a  Channel  Title,  “chtitle”.  These  assignments 
will  be  utilized  later  for  labeling  the  generated  jpg  figures. 
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%  Looping/ Processing  of  Data  Channels 
for  ch  =  1: channels 

if  channels==46 


f scaler 
psdunits 
raounits 
chtitle 
el seif  ch==2 
f scaler 
psduni ts 
raouni ts 
chtitle 


fsl; 
pul; 
rul ; 

'Wave  Ht  Bow' ; 

fsl; 

pul; 

rul; 

' Sonix  Sonic'; 


else  if  ch==46 
f scaler  = 

psdunits  = 
raouni ts  = 
chtitle  = 
else 
end 

elseif  channels==38 

if  ch-=l 

fscaler  = 
psduni ts  = 
raounits  = 
chtitle  = 
elseif  ch==2 
fscaler  = 
psduni ts  = 
raounits  = 
chtitle  = 


fs4; 
pu4; 
ru  4; 

’Lghtrl-PBow  TvAcc'; 


fsl; 

pul; 

rul; 

'Wave  Ht  Bow' ; 

fsl; 

pul; 

rul; 

'Sonix  Sonic'; 


elseif  ch==38 

fscaler  =  fs4; 
psduni  ts  =  pu4  ; 
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raounits  =  ru4; 
chtitle  =  'Lghtrl-PBowTvAcc'; 
else 
end 

elseif  channels==24 


f scaler 
psduni ts 
raounits 
chtitle 
elseif  ch==2 
f scaler 
psdunits 
raounits 
chtitle 


fsl ; 
pul; 
rul; 

'Wave  Ht  Bow' ; 

fsl; 

pul; 

rul; 

' Sonix  Sonic'; 


elseif  ch==24 
fscaler  - 
psduni ts  = 
raounits  = 
chtitle  = 
else 
end 


fs4; 

pu4; 

ru4; 

' BoomTip-Vert  Acc '  ; 


end 


8.  Base-lining  the  Raw,  Time-based  Data 

The  actual  processing  of  data  begins  here  by  base-lining  or  removing  the  mean 
from  the  raw,  time-based  data.  Otherwise,  “Failure  to  remove  large  sample  means  or 
other  trends  in  the  data  may  result  in  distorted  or  biased  spectral  estimates.” 
(Marple,  1 987, pp.  1 32) 
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%  Base-lining  the  Raw,  Time-based  Data - 

x  =  data ( : , ch) ; 

x  =  x-mean  (x)  ; 

9.  FFT  Computation 

The  MATLAB  command  “fft”  returns  the  Fourier  transform  of  the  base-lined, 
column  data  matrix.  If  the  matrix  length  is  a  power  of  two,  the  “fft”  command  employs  a 
high-speed  fast  Fourier  transform  algorithm.  But,  this  is  not  always  the  case  since  the 
number  of  data  points  per  channel  varies  from  run  to  run  In  this  situation,  an  alternate 
mixed-radix  algorithm  finds  the  prime  factors  of  the  column  matrix  length  then  computes 
the  discrete  Fourier  transforms  of  the  shorter  sequences.  Either  way,  the  Fourier 
transform  of  the  base-lined  data  is  produced. 


%  FFT  Computation 
Y  =  fft(x,n); 


10.  Model-Scale  PSD  Computation 

The  Model-Scale  PSD  is  computed  by  multiplying  the  elements  of  the  Fourier 
Transform  by  their  complex  conjugate  and  dividing  by  the  total  number  of  elements. 


%  Model-Scale  PSD  Computation 
PSD  =  Y.  *conj  (Y) /n; 


21 


11.  Movir  ^-window  Averaging  of  Model-Scale  PSDs 


This  is  where  the  RAO  program  method  differs  significantly  from  the  tools  in 
MATLAB  previously  identified  by  offering  a  much  finer  control  over  the  resolution  vs. 
display  accuracy  issue.  An  envelope  initially  starts  with  half  the  “ptavg”  window  plus 
one  PSD  points,  averages  them  together  and  assigns  that  average  as  the  first,  smoothed 
PSD  point.  The  envelope  then  accepts  the  next  sequential  PSD  point,  averages  the  now 
larger  window  and  assigns  that  average  as  the  second,  smoothed  PSD  point.  The  process 
continues  until  the  envelope  reaches  it’s  maximum  size  of  “ptavg”  plus  one  where  the 
latest  average  is  assigned  as  the  smoothed  PSD  point  for  the  middle  of  the  envelope. 
From  here  on  out,  the  envelope  accepts  the  next  sequential  PSD  point  while  dropping  the 
first  one,  thus  maintaining  a  fixed-length,  moving  window.  Refer  to  the  following 
example: 
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Example:  Let  “ptavg”  =  4 


[ . Freshly.  Produced  PSD  Data  Points . 1 

HI  #2  #3  #4  #5  #6  #7  #8  #9  #M 

[  1st  avg  PSD  ]  Window  starts  at  half  “ptavg  ”  plus  one  (i.e.  2  +  1) 

[  2nd  avg  PSD  ]  Window  increments  in  size  by  one 

[  3rd  avg  PSD  ]  Max  capacity  of  ptavg  +  1  (i.e.  4  +  1) 

[  4th  avg  PSD  ]  Window  begins  shift  to  the  right 

[  5th  avg  PSD  ]  Window  continues  to  shift 

[  6th  avg  PSD  ]  Window  still  shifting 

[  .  7 

At  this  point  in  the  RAO  program,  the  smoothed  PSDs  are  designated  as  a  column 
of  data  in  a  composite  PSD  matrix  with  the  column  designation  correlating  directly  to  the 
data  channel  being  processed.  Thus,  the  first  channel  processed  becomes  the  first  column 
of  PSD  data  in  the  composite  PSD  matrix.  As  for  the  chosen  window  length  of  50  plus 
one,  tests  of  trial  runs  using  incrementally  larger  point  average  envelopes  yielded  a 
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compromise  of  “ptavg”  =  50  which  offered  the  best  display  resolution  while  maintaining 
accuracy  in  comparison  to  the  spreadsheet  results  provided  (see  figs.  1-3) 


%  Moving-window  Averaging  of  Model-Scale  PSDs  for 
Smoothing 

halfptavg  =  ptavg/2 ; 
sum  =  0; 
for  k  =  1: ptavg 

sum  =  sum  +  PSD(k,  1) ; 

if  k  >  halfptavg 

psd (k-half ptavg, ch)  =  sum/k; 
else 
end 

end 

pts  =  window  +  halfptavg; 

ptavgplusone  =  ptavg  +  1; 

j  =  1; 

for  k  =  ptavgplusone : pts 

sum  =  sum  +  PSD  (ptavg+j ,  1)  -  PSD(j,l); 

psd ( k-half ptavg , ch)  =  sum/ptavg; 

j 

end 

meansqr 
psdarea 
psdscaler 
psd ( : , ch) 


=  j  +  1; 


=  mean  (x.  *2)  ; 

=  trap  (momega,psd( : ,  ch)  )  ; 
=  meansqr /psdarea; 

=  psd  ( : ,  ch)  *psdscaler 
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[Spreadsheet  RAO 
•Results  V 


Figure  3. 


12.  Computation  of  Model-Scale  RAOs 

As  previously  addressed  in  the  theory,  RAOs  are  obtained  by  taking  the 
individual  square  roots  of  the  desired  PSDs  divided  by  the  seaway  or  wave  height  PSDs. 
With  the  first  data  channel  for  each  configuration  allocated  as  the  reference  wave  height, 
the  first  channel  of  RAOs  processed  will  naturally  be  at  unity,  a  constant  dimensionless 
value  of  one. 


%  Computation 
X 

mrao ( : , ch) 


of  Model-Scale  RAOs - 

=  psd  ( : ,  ch)  .  /psd  (:,  1)  ; 
=  sqrt  (X)  ; 
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13.  Computation  of  Full-Scale  RAOs 


Computation  of  the  Full-Scale  RAOs  is  simply  a  process  of  multiplying  the 
Model-Scale  RAOs  by  the  Full-Scale  Scaler.  This  previously  selected  scaler  is  based 
upon  the  nature  of  the  data  sampled. 


%  Computation  of  Full-Scale  RAOs - 

frao(:,ch)  =  mrao  (:,  ch)  *f scaler; 


14.  Plot  Model-Scale  PSDs  and  Save  Figure  in  jpg  Format 

A  plot  of  the  smoothed,  Model-Scale  PSDs  vs.  their  associated  Model-Scale 
omegas  is  generated  (see  fig.  4)  then  saved  as  a  figure  in  jpg  format  in  a  folder  dedicated 
to  the  particular  run  and  titled  accordingly  with  the  run  number,  channel  number  and  a 
“p”  for  Model-Scale  PSDs  label. 


%  Plot  Model-Scale  PSDs  and  Save  Figure  in  jpg  Forma  t- 
clf 

plot  (momega  (1,  :)  ,  psd  ( : ,  ch)  ) 
grid 

sch  =  num2str(ch); 

title ([ 'Model-Scale  PSDs  for  DTMB  Run  'srun,  '  / 

Channel  ',sch,':  'chtitle]) 
xlahel ( 'Omega  (rad/ s)  ' ) 
yl abel ( [ ' PSD  'psduni ts] ) 

saveas  (gcf,  [swrite,  folder ,  '  \ ' ,  filename,  'C ,sch,  'p'J, 

'jpg’); 
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Figure  4. 


15.  Plot  Model-Scale  RAOs  and  Save  Figure  in  jpg  Format 

A  plot  of  the  Model-Scale  RAOs  vs.  their  associated  Model-Scale  omegas  is 
generated  (see  fig.  5)  then  saved  as  a  figure  in  jpg  format  in  a  folder  dedicated  to  the 
particular  DTMB  run  and  titled  accordingly  with  the  run  number,  channel  number  and  an 
“m”  for  Model-Scale  RAOs  label. 


%  Plot  Model-Scale  RAOs  and  Save  Figure  in  jpg  Format  - 
clf 

plot  (momega  (1,  : ) , mrao  ( : ,  ch) ) 
grid 

title  ([  'Model-Scale  RAOs  for  DTMB  Run  'srun,  '  / 

Channel  ' , sch,  ' :  'chtitle]) 
xlabel ( 'Omega  (rad/s) ' ) 
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ylabel ( [ 'RAO  'raouni ts] ) 

saveas  (gcf,  [swrite,  folder,  '  \ filename,  'C',sch,  'm'], 

'jpg') ; 


Figure  5. 


16.  Plot  Full-Scale  RAOs  and  Save  Figure  in  jpg  Format 

A  plot  of  the  Full-Scale  RAOs  vs.  their  associated  Full-Scale  omegas  is  generated 
(see  fig.  6)  then  saved  as  a  figure  in  jpg  format  in  a  folder  dedicated  to  the  particular  run 
and  titled  accordingly  with  the  run  number,  channel  number  and  an  “f’  for  Full-Scale 
RAOs  label. 
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%  Plot  Full-Scale  RAOs  and  Save  Figure  in  jpg  Format— 
elf 

plo t  (  f  omega  ( 1 ,  : ) ,  frao  ( : ,  ch)  ) 
grid 

title ([' Full-Scale  RAOs  for  DTMB  Run  'srun,  '  / 

Channel  ',sch,  '  chtitle] ) 
xlabel (' Omega  (rad/s)') 
ylabel ( [ 'RAO  ' raounits] ) 

saveas  (gcf,  [swrite,  folder ,  '  \  filename,  'C'  ,sch,  'f' ], 
j PS ' )  ; 


Processing  of  the  first  data  channel  is  now  complete.  The  RAO  program  then 
loops  back  to  event  7  listed  above,  Looping/Processing  of  Data  Channels,  until  all  data 
channels  for  the  specific  run,  as  determined  by  the  data  configuration,  have  been 
processed  to  generate  the  Model-Scale  PSDs,  Model-Scale  RAOs  and  Full-Scale  RAOs. 


17.  Save  Omega  Matrices  as  Last  Column  of  Composite  PSD  and  RAO 
Matrices 

For  the  purpose  of  convenient  analysis,  the  single-column  omega  matrices  are 
tacked  on  as  the  final  columns  of  their  respective  PSD  or  RAO  composite  matrices.  This 
basically  eliminates  the  necessity  of  the  jpg  formatted  figures  by  allowing  direct 
comparison  between  any  two  columns  of  a  composite  matrix  simply  by  using  a  plotting 
function.  Thus,  PSD  vs.  PSD,  PSD  vs.  omega,  RAO  vs.  RAO  and  RAO  vs.  omega  plots 
can  readily  be  produced 


%  Save  Omega  Matrices  as  Last  Column  of  Composite  PSD 
and  RAO  Matrices 
ch  =  ch  +  1; 

psd ( : , ch)  =  momega  ( 1 , : ) ' ; 
mrao  ( : ,  ch)  =  momega  (1,  :)  '  ; 
frao ( : , ch)  =  f omega (1,  : ) ' ; 


18.  Save  Composite  PSD  and  RAO  Matrices 

This  final  event  in  the  sequence  saves  the  composite  PSD  and  RAO  matrices 
along  with  the  jpg  formatted  figures  previously  saved  in  the  folder  dedicated  to  the 
particular  DTMB  run.  It’s  now  time  to  process  another  DTMB  run. 


%  Save  Composite  PSD  and  RAO  Matrices - 

save(  [swrite,  folder,  '\\  filename,  'psd.txt'],  'psd', 
'-ascii' , '-tabs'); 

save  ( [swrite,  folder ,  '  \ ',  filename,  'mrao.txt'],  'mrao', 
' -ascii' , '-tabs'); 

save ( [swrite, folder,  ' \ ', filename, 'frao.txt'], 'frao', 
' -ascii ' , ' -tabs ' ) ; 
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B.  VALIDATION 


Validation  of  the  RAO  Program  was  performed  by  utilizing  a  MATLAB 
comparison  program,  “psdandraocomparisons.m”,  written  to  compare  the  various 
methods  of  resultant  PSDs;  Welch,  Periodogram,  Multi-taper,  Spreadsheet  and  RAO 
Program.  The  comparison  program  is  comprised  of  a  sequence  of  events  with  each 
explained  individually  as  applied  to  run  67: 


1.  Load  Applicable  Matrices 

Once  the  comparison  program  has  been  initiated,  the  following  matrices  are 

loaded: 


Matrix: _ 

R067 

Run67.txt 

R067psd.txt 

R067mrao.txt 

R067frao.txt 


Definition: _ 

Raw,  Time-based  Data  (All  46  Channels) 
Spreadsheet  Model-Scale  Omegas,  PSDs  and  RAOs 
RAO  Program  Model-Scale  PSDs  and  Omegas 
RAO  Program  Model-Scale  RAOs  and  Omegas 
RAO  Program  Full-Scale  RAOs  and  Omegas 
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%  Load  Applicable  Matrices 

load  R067; 

load  Run67.txt; 

load  R067psd.  txt; 

load  R067mrao.txt; 

load  R067frao.txt; 


2.  Base-line  the  Raw  Time-Based  Data 

Prior  to  initiating  the  MATLAB  methods  of  Welch,  Periodogram  and  Multi-taper, 
the  raw,  time-based  data  for  both  Channel  1:  Wave  Ht  and  Channel  3:  TACS  Roll  are 
base-lined  for  continuity  just  as  it  is  done  in  the  RAO  Program. 


%  Base-line  the  Raw  Time-based  Data  from  R06 7 


X 

=  data(:,l); 

% 

Wave 

Ht  Data 

X 

=  x  -  mean  (x)  ; 

z 

=  data  ( : , 3) ; 

% 

TACS 

Roll  Data 

z 

=  z  -  mean(z) ; 

3.  Compute  Wave  Ht  PSDs  using  MATLAB  Methods 

The  Welch,  Periodogram  and  Multi-taper  methods  are  used  to  compute  their 

Wave  Ht  PSDs  and  associated  Frequencies.  The  Frequencies  will  be  converted  to 

Omegas  later.  For  efficiency  of  computation,  a  maximum  of  8096  FFTs  is  designated. 

This  happens  to  be  the  largest  power  of  two  that  is  less  that  the  total  number  of  data 

points  per  channel  which  range  from  approximately  14,000  to  16,000.  Otherwise, 

MATLAB  uses  a  slower  approach  of  seeking  out  the  primes  of  the  total  number  of  data 

points  and  computing  the  FFTs  according  to  these  segments. 
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%  Computation  of 
[hwel ch , fwel ch ] 
[hperio,  fperio] 
[hmul  t i ,  fmul  ti] 
meansqrx 
wel chareah 
perioareah 
multiareah 
hwelchscaled 
hperioscaled 
hmultiscaled 


Wave  Ht  PSDs  and  Associated  Freqs- 
=  pwelch  (x,  [],[],  8096 ,  32 . 2 )  ; 

=  periodogram  (x,  [] ,  8096, 32 . 2)  ; 

=  pmtm  (x,  [ ]  ,  8096,  32 . 2) 

=  mean  (x.  ^2)  ; 

=  trap  ( wwel  ch ,  hwel  ch )  ; 

=  trap  ( wperi o ,  hperi o)  ; 

=  trap  (wmulti,hmulti) ; 

=  hwelch  *meansqrx/wel  chareah  ; 

=  hperi  o  *meansqrx/peri  oareah  ; 

=  hmul  ti  *meansqrx/mul  tiareah; 


4.  Compute  TACS  Roll  PSDs  using  MATLAB  Methods 

The  same  process  described  above  is  now  performed  for  the  TACS  Roll  PSDs. 


%  Computation  of 
[  rwel  ch ,  fwel  ch  ] 

[ rperi o ,  fperi o] 
[rmulti, f multi] 
meansqrz 
welcharear 
perioarear 
multiarear 
jrwelchscaled 
rperioscaled 
multiscaled 


TACS  Roll  PSDs  and  Associated  Freqs - 
=  pwelch  (z,  [],[],  8096,  32 .2)  ; 

=  periodogram  ( z,  [],  8096,32 .2)  ; 

=  pmtm(z,  [ ]  ,  8096,  32 .2)  ; 

=  mean  (z.  *2)  ; 

=  trap  ( wwel ch , rwel ch ) ; 

=  trap ( wperi o , rperi o) ; 

=  trap  (wmulti,  multi)  ; 

=  rwelch*meansqrz /welcharear  ; 

=  rperio  *meansqrz /perioarear  ; 

=  rmul ti *meansqrz /mul tiarear; 


5.  Convert  MATLAB  Method  Frequencies  to  Omegas 

Conversion  of  frequency  to  omega  is  done  simply  by  multiplying  the  frequency 
matrices  by  2 it. 
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%  Convert  Frequency  to  Omega - 

wwel ch  =2 *pi *fwel ch ; 

wperi  o  =  2  *pi  *fperi  o  ; 

wmulti  =  2*pi  * f multi; 


6.  Compute  Areas  under  the  PSD  Curves 

As  an  integrity  check  along  the  way,  the  mean  squares  of  the  base-lined  data  are 
compared  to  the  areas  under  the  PSD  curves.  At  this  point,  both  the  Spreadsheet  and 
RAO  Program  results  are  integrated  into  the  comparison  program  to  ensure  continuity 
between  the  various  methods  with  comparative  results  as  follows: 


a.  Wave  Ht  Comparison 

Method: _ Result; 


Mean  Square  of  Base-lined  Wave  Ht  Data 

0.0124 

Welch  Area  (rad/s) 

0.0124 

Periodogram  Area  (rad/s) 

0.0124 

Multi-taper  Area  (rad/s) 

0.0124 

Spreadsheet  Area  (rad/s) 

0.0126 

RAO  Program  Area  (rad/s) 

0.0124 
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From  the  results  listed  above  for  run  67,  the  Welch,  Periodogram,  Multi-taper  and 
RAO  Program  methods  agree  to  the  nearest  l/1000th  when  rounded  to  5  decimal 
positions.  In  addition,  each  is  within  1.6%  of  the  Spreadsheet  results. 


b.  TACS  Roll  Comparison 

Method: 

Result: 

Mean  Square  of  Base-lined  TACS  Roll  Data 

0.0061 

Welch  Area  (deg/s) 

0.0061 

Periodogram  Area  (deg/s) 

0.0061 

Multi-taper  Area  (deg/s) 

0.0061 

Spreadsheet  Area  (deg/s) 

0.0059 

RAO  Program  Area  (deg/s) 

0.0061 

From  the  results  listed  above  for  run  67,  the  Welch,  Periodogram,  Multi-taper  and 
RAO  Program  methods  agree  to  the  nearest  1/1 000th  when  rounded  to  5  decimal 
positions.  In  addition,  each  is  within  3.3%  of  the  Spreadsheet  results. 
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%  Compute  Areas 
welchareah2 
perioareah2 
multiareah2 


under  the  Wave  Ht  PSD  Curves -- 
=  trap  (wwelch,  hwelchscaled)  ; 
=  trap  (wperio, hperioscaled)  ; 
=  trap  (wmulti ,  hmultiscaled)  ; 


spreadareah  =  trap  (run67  ( : ,  2)  ,  run67  ( : ,  3)  )  ; 

raoprogareah  =  trap  (R067psd ( : ,  47) , R067psd (:,!)); 

%  Compare  Areas  under  the  Wave  Ht  PSD  Curves  against 
meansqrs 

wavehtpsdareacomparison  =  [meansgrx  welchareah2 

perioareah2  multiareah2  spreadareah  raoprogareah] 


%  Compute  Areas  under  the  TACS  Roll  PSD  Curves 
welcharear2  =  trap  (wwelch,  rwelchscaled)  ; 

perioarear2  =  trap  (wperio,  rperioscaled) ; 

multiarear2  =  trap (wmulti, rmultiscaled); 


spreadarear  =  trap (run67 (:  ,2) , run67 (:,  4) ) ; 
raoprogarear  =  trap  (R067psd  (:  ,47)  ,  R067psd  (:  ,3)  )  ; 

%  Compare  Areas  under  the  TACS  Roll  PSD  Curves  against 
meqnsqrs 

rollpsdareacomparison  =  [meansqrz  welcharear2 

perioarear2  multiarear2  spreadarear  raoprogarear] 


7.  Graphical  Comparison  of  Model-Scale  PSDs 

Plots  of  Model-Scale  PSDs  from  all  five  methods  for  both  the  Wave  Ht  and 
TACS  Roll  are  generated  for  visual  comparison.  From  these  comparison  plots,  it  was 
determined  that  previously  discussed  “ptavg”  selection  of  50  was  a  very  good 
compromise  between  display  resolution  and  accuracy  when  compared  to  the  other  four 
methods  of  generating  PSDs.  (see  figs.  7-10) 
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%  Graphical  Model-Scale  Wave  Ht  PSDs  Comparison - 

All  5  Methods 
figured) 

plot (wwelch, hwelchscaled)  %  Welch  Method 

hold  on 

plot (wperio, hperioscaled)  %  Perio  Method 

plot (wmulti,hmultiscaled)  %  Multi  Method 

plot (run67 ( : , 2 ) , run67 ( : , 3 ) )  %  Spreadsheet 

plot  (R067psd  (:,47),  R067psd  (:,!));  %  Rao  Program 

title  ( 'R067CH01  Model-Scale  Wave  Ht  PSDs  Comparison' ) 
xlabel (' Omega  (rad/s)') 
ylahel ( 'PSD  (dimesionless) ' ) 
axis ( [0  6  0  .015]) 
grid 

%  Graphical  Model-Scale  Wave  Ht  PSDs  Comparison - 

Spreadsheet  vs.  RAO  Program  only 
figure (2) 

plot (run67 ( : , 2) , run67 ( : , 3) )  %  Spreadsheet 

hold  on 

plot (R067psd ( : , 47 ) , R067psd ( : , 1 ) )  %  Rao  Program 

title ( 'R067CH01  Model-Scale  Wave  Ht  PSDs  Comparison' ) 

xlabel ( 'Omega  (rad/s) ' ) 

ylabel ( 'PSD  (dimesionless ) ' ) 

axis ( [0  6  0  .0065]) 

grid 

%  Graphical  Model-Scale  TACS  Roll  PSDs  Comparison - 

All  5  Methods 
figured) 

plot (wwelch, rwelchscaled)  %  Welch  Method 

hold  on 

plot (wperio,  rperioscaled)  %  Perio  Method 

plot (wmulti,  rmultiscaled)  %  Multi  Method 

plot (run67 ( : , 2 ) , run67 ( : , 4) )  %  Spreadsheet 

plot (R067psd( : , 47 ) ,R067psd( : , 3) ) ;  %  Rao  Program 

title ( 'R067CH03  Model-Scale  TACS  Roll  PSDs 
Comparison ' ) 
xlabel ( 'Omega  (rad/s)') 
ylabel ( 'PSD  (deg)  ') 
axis ( [0  60  .02]) 
grid 
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%  Spreadsheet 


%  Graphical  TACS  Roll  PSDs  Comparison- 
Spreadsheet  vs.  RAO  Program  only 
figure  (4) 

plot  (run67  ( :  ,2)  ,  run67  ( : ,  4)  ) 
hold  on 

plot  (R067psd  ( :  /  47)  ,R067psd  (:  ,3)  )  %  Rao  Program 

title ( 'R067CH03  Model-Scale  TACS  Roll  PSDs 
Comparison ' ) 
xlabel ( 'Omega  (rad/s) ' ) 
y label  ( 'PSD  (deg)') 
axis ([060  . 009] ) 
grid 


Figure  7. 
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Figure  10. 


8.  Compute  Model-Scale  RAOs 

Just  as  it  is  performed  in  the  RAO  Program,  Model-Scale  RAOs  are  computed  for 
the  Welch,  Periodogram  and  Multi-taper  methods  by  taking  the  square  of  the  scaled 
TACS  Roll  PSDs  divided  by  the  scaled  Wave  Ht  PSDs. 


%  Compute 
Rl 

welchmrao 

R2 

periomrao 

R3 

mul timrao 


Model-Scale  TACS  Roll  RAOs - 

=  rwelchscaled. /hwelchscaled; 
=  sqrt  (Rl )  ; 

=  rperioscaled. /hperioscaled; 
=  sqrt  (R2)  ; 

=  rmultiscaled. /hmultiscaled; 
=  sqrt(R3); 
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9.  Graphical  Comparison  of  Model-Scale  RAOs 

Plots  of  Model-Scale  TACS  Roll  RAOs  from  all  five  methods  are  generated  for 
visual  comparison  and  they  provide  very  satisfactory  results,  (see  figs.  11-12) 


Figure  11. 


Figure  12. 


10.  Scale  to  Full-Scale  RAOs  and  Omegas 

In  order  to  compare  Full-Scale  RAOS  and  Omegas,  the  Model-Scale  RAOs  and 
Omegas  are  scaled  accordingly  just  as  they  are  done  in  the  RAO  Program. 


%  Scale  to 
welch frao 
periofrao 
multifrao 


Full-Scale  RAOs  and  Omegas - 

=  welchmrao*12/24 . 175; 
=  periomrao* 12/24 . 175 ; 
=  multimrao*12/24 . 175 ; 


wwelchfull 

wperiofull 

wmultifull 


wwelch/sqrt (24.1 75) ; 
wperio/sqrt  (24.175)  ; 
wmulti/sqrt (24.1 75) ; 
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11.  Graphical  Comparison  of  Full-Scale  RAOs 

Finally,  the  Full-Scale  RAOs  and  Omegas  for  the  Welch,  Periodogram,  Multi¬ 
taper  and  RAO  Program  methods  are  plotted  for  visual  comparison,  (see  figs.  13-14) 


%  Graphical  Full-Scale  TACS  Roll  RAOs  Comparison - 

4  Methods  Not  including  Spreadsheet 
figured) 

plot (wwelchfull ,welchfrao)  %  Welch 

hold  on 

plot  (vjperio full, per iofrao)  %  Periodogram 

plot (wmulti full, multi frao)  %  Multi-taper 

plot  (R067frao(: ,  47)  ,R067frao(:  ,3))  %  RAO  Program 

title  (' R067CH03  Full-Scale  TACS  Roll  RAOs  Comparison' ) 
xlabel  (' Omega  (rad/s)') 
ylabel ( 'RAO  ( deg/ in) ' ) 
axis ( [0  1.4  0  3]) 
grid 

%  Graphical  Full-Scale  TACS  Roll  RAOs  Presentation - 

RAO  Program  only- 
figure  (8) 

plot  (RO  67  frao  (:,47)  ,  R067frao  (:,3)  )  %  Rao  Program 

title ( 'R067CH03  Full-Scale  TACS  Roll  RAOs 
Presentation ' ) 
xlabel ( 'Omega  (rad/s) ' ) 
ylabel ('RAO  (deg/in)') 
axis ( [0  1.4  0  1.4]) 
grid 
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IV.  CONCLUSIONS  &  RECOMMENDATIONS 


A.  CONCLUSIONS 

Due  to  the  validation  comparisons  performed  on  run  67  above,  the  RAO  Program 
performed  quite  admirably  and  accurately.  The  RAO  Programs  ability  to  accurately 
produce  Model-Scale  PSDs,  Model-Scale  RAOs  and  ultimately  Full-Scale  RAOs  for 
each  channel,  along  with  their  associated  Omegas,  attests  to  the  programs  performance. 
In  addition,  a  short  looping  program,  “raoloop.m”,  enabled  up  to  seven  data  runs  to  be 
evaluated  sequentially. 


%  Short  Program  Title:  "raoloop.m" - 

a  =  [19  21  23  25  27  29  31] ;  %  DTMB  Run  Numbers 

for  z  =  1:7 

rao  (a(l,  z)  )  ; 

end 


The  seven  run  looping  limit  was  imposed  due  to  local  hardware  configuration 
limitations.  With  each  data  run  producing  approximately  12.5  megabytes  of  stored 

l 

figures  and  matrices,  a  100  megabyte  zip  disc  could  safely  store  up  to  seven  data  runs. 
With  a  read/write  CD  drive,  the  looping  capability  would  have  been  much  greater. 


Upon  completion  of  processing  the  Configuration  I  runs,  the  composite  matrices 
(Model-Scale  PSDs,  Model-Scale  RAOs  and  Full-Scale  RAOs)  for  each  run  were 
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assembled  all  on  one  compact  disc  as  a  convenient  reference  from  which  any  two 
columns  of  results  may  be  plotted  against  each  other. 


B.  RECOMMENDATIONS 

An  opportunity  for  a  logical  continuation  of  this  thesis  will  soon  be  available. 
Trial  runs  using  the  actual,  Full-Scale  ship,  the  USNS  Grand  Canyon  State  (ex-SS 
President  Polk)  (T-ACS  3),  were  conducted  September  9th- 1 6th  off  the  coast  of  Southern 
California  near  Camp  Pendleton.  Data  recorded  from  the  Full-Scale  trial  runs  could  be 
processed  in  a  manner  similar  to  the  RAO  Program  and  used  as  a  comparison  tool  to 
either  disprove  or  further  prove  the  validity  of  results  from  the  RAO  Program.  If  RAO 
results  from  both  the  at-sea  Full-Scale  trial  runs  and  the  RAO  Program  model-basin 
evaluations  are  in  agreement,  a  multitude  of  combinations  can  then  be  evaluated  to 
predict  ship  motions  in  an  irregular  sea  by  using  the  principle  of  linear  superposition. 
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APPENDIX  A:  TEST  CONFIGURATION  I 
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TASCA5  1395,397  1169 


TASCA5  13  97,399  1170 


400,402  1171 


402,404  1172 


404,406  1173 


406,408  1174 


409,411  1175 


411,413  1176 


ASCA5  1413,415  1177 


ASCA6  1416,418  1178 


TASCA6  1418,420  1179 


TASCA6  1420,422  1180 


TASCA6  1422,424  1181 


TASCA6  1424,426  1182 


429,431  1183 


270  13s  1270  136.5  1  129.6 


443,445  118811  1180 


445,447  1189  ]l  1180  14s  270  136.5  12  154.5 


180  13  s  1270  136.5  1  129.6 


24  file  format. 


270  12  1270  136.5  112  154.5 


270  13  1270  136.5  1  154.5 


270  13s  1270  136.5  I  154.5 


270  4  1270  136.5  112  154.5 


270  12  1270  136.5  1  154.5 


270  13  1270  136.5  1  154.5 


270  13  s  1270  136.5  1  154.5 


270  |2  10  136.5  1  154.5 


270  13  10  13  6.5  I  154.5 


270  |3s  10  136.5  15.5  154.5 


270  l3s  |0  136.5  1  154.5 


270  13  10  136.5  112  154.5 


180  12  1270  13  6.5  I  154.5 


TASCA7  1431,433  184  I  1180  3  1270  136.5  |12  154.5 


TASCA7  1433 , 435  Il85il  1270  3  270  136.5  12  154.5 


435,437  1186 1!  1180  J3s  1270  136.5  112  154.5 


437,439  §187 II  180  14  1270  J36.5  112  154.5 


270  136.5  1  154.5 


24  File  format 
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490 

15- 

Aug- 

97 

TASCA7 

489 , 491 

211 

I 

270 

3 

315 

29.25 

12 

25 

492 

15- 

Aug- 

97 

flU 

491,493 

212 

I 

180 

3s 

315 

29.25 

12 

25  j 

494 

15- 

Aug- 

97 

SSI 

493,495 

213 

I 

270 

3s 

315  | 

29.25  | 

12 

25 

496 

a 

mm 

SSSi 

495,497 

214) 

I 

180  : 

4 

315  | 

29.25 

12 

25 

■■ 

498 

15- 

Aug- 

97 

51 

497,499 

215| 

i  S 

270  | 

4 

315  | 

29.25 

12 

25 

H 

500 

15- 

Aug- 

97 

TASCA7 

499,501 

216 

_ j 

1 

180 

4s 

315 

29.25 

12 

25 

■H 

502 

15-  | 

Aug- 

97 

ssss 

501 

217 

i 

270  | 

4s  ! 

315  i 

29.25 

12 

25 
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APPENDIX  B.  TEST  CONFIGURATION  H 


503 , 505  |218ili  1180  |3  1315  15 


TACSA7  *505 , 507  1219 1  1270  I  1315 


TACSA7  1507,509  1220 


516,518  1224 


ASCA8  1525 


270  |3s  i  15 


509,511  12 2 1  jl i  1270  f  1315  15 


511,513  1222  Jli  1270  14s  1315  15 


90  I  1270  I 


90  J3s  1270  15 


518,520  225  jli  190  14  1270  |5 


520,522  226) Ii  j  90  4s  1 270  5 


522,524  223 1  Ii  90  3  1 270  5 


CD  Ii  90  4 


CD  Ii  (90  4 


CDili  190  4 


CD  Jli  190  14 


12 

25 

12 

25 

12 

25 

12 

25 

12 

25 

12 

25 

12 

25 

12 

25 

12 

25 

12 

25  | 
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APPENDIX  C:  TEST  CONFIGURATION  III 


9 

1 

1 

■ 

a 

ai 

■ 

■■ 

■ 

■ 

576 

17- 

Aug- 

97 

TASCA9 

575 , 577 

244 

iii 

0 

3s 

90 

45.3 

13 

60  j 

578 

TASCA9 

577,579 

245 

... 

m 

90 

3s 

90 

45.3 

13 

60  j 

580 

H 

■Sm 

579,582 

246 1 

. 

m 

0 

4 

90 

45.3 

13 

60 

581 

s 

TASCA9 

579,582 

247! 

■ 

0 

4  + 

90 

45.3 

13 

60 

583 

i 

B 

582 

248 1 

m 

90 

3 

90 

45.3 

13 

60 

240secs  Static  + 
Load  Control 

Demo 

584 

9 

mm 

582  | 

249 

iii 

0 

2 

90 

0 

0 

60 

240secs  Static  + 
Load  Control 

Demo 

585 

a 

mm 

582 

250 

iii  j 

90 

2 

90 

0  | 

0 

60  j 

586 

a 

IBB 

TASCA9 

582 

251 

90 

3 

90  j 

0 

0 

60  1 
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APPENDIX  D:  MASTER  CHANNEL  LIST:  CONFIGURATION  I 


DTMB  Runs  1-503  and  NRL  Runs  1-217 


Channel  # 

Channel  Title 

01 

Wave  Ht  Bow 

02 

Sonix  Sonic 

03 

TACS  Roll 

04' 

TACS  Pitch 

05 

TACS  Roll  Rt 

06 

TACS  Pitch  Rt 

07 

TACS  Yaw  Rt 

08 

TACS  CG  ZAcc 

09 

TACS  CG  YAcc 

10 

TACS  CG  XAcc 

11 

TACS  CT  ZAcc 

12 

TACS  CT  YAcc 

13 

Cntr  Roll 

14 

Cntr  Pitch 

15 

Cntr  CG  ZAcc 

16 

Cntr  CG  YAcc 

17 

Cntr  CG  XAcc 

18 

TACS  Cntr  Rel  X 

19 

TACS  Cntr  Rel  Y 

20 

TACS  Cntr  Rel  Z 

21 

Lhtr  Roll 

22 

Lhtr  Pitch 

23 

Lhtr  Roll  Rt 

24 

Lhtr  PitchRt 

25 

Lhtr  CG  ZAcc 

26 

Lhtr  CG  YAcc 

27 

Lhtr  CG  XAcc 

28 

TACS  Lhtr  Rel  X 

29 

TACS  Lhtr  Rel  Y 

30 

TACS  Lhtr  Rel  Z 

31 

Load  Mo  Orthognl 

32 

Load  Mo  Parallel 

33 

BoomTip-HorizAcc 

34 

BoomTip-Vert  Acc 

35 

Lghtr4-Strn  VAcc 

36 

Lghtr4-Bow  VAcc 

37 

Lghtr3-Bow  VAcc 

Units 

Final  File  Title 

inch 

r###ch01 

inch 

r###ch02 

degree 

r###ch03 

degree 

r###ch04 

deg/sec 

r###ch05 

deg/sec 

r###ch06 

deg/sec 

r###ch07 

& 

r###ch08 

G 

r###ch09 

G 

r###chl0 

G 

r###chl  1 

G 

r###chl2 

degree 

r###chl3 

degree 

r###chl4 

G 

r###chl5 

G 

r###chl6 

G 

r###chl7 

inch 

r###chl8 

inch 

r###chl9 

inch 

r###ch20 

degree 

r###ch21 

degree 

r###ch22 

deg/sec 

r###ch23 

deg/sec 

r###ch24 

G 

r###ch25 

G 

r###ch26 

G 

r###ch27 

inch 

r###ch28 

inch 

r###ch29 

inch 

r###ch30 

inch 

r###ch31 

inch 

r###ch32 

G 

r###ch33 

G 

r###ch34 

G 

r###ch35 

G 

r###ch36 

G 

r###ch37 
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38 

Lghtr2-Bow  VAcc 

G 

r###ch38 

39 

Lghtrl-Bow  VAcc 

G 

r###ch39 

40 

Lghtr2-SMdShpVAc 

G 

r###ch40 

41 

Lghtr2-PMdShpVAc 

G 

r###ch41 

42 

Lghtr2-CLMdSpVAc 

G 

r###ch42 

43 

Lghtrl  -CLMdSpVAc 

G 

r###ch43 

44 

Lghtr4-StrnLngAc 

G 

r###ch44 

45 

Lghtr4-PStrnTvAc 

G 

r###ch45 

46 

Lghtrl  -PBowTvAcc 

G 

r###ch46 
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APPENDIX  E.  MASTER  CHANNEL  LIST:  CONFIGURATION  n 


DTMB  Runs  504-533  and  NRL  Runs  21 8-226 


Channel  # 

Channel  Title 

Units 

Final  File  Title 

01 

Wave  Ht  Bow 

inch 

r###ch01 

02 

Sonix  Sonic 

inch 

r###ch02 

03 

TACS  Roll 

degree 

r###ch03 

04 

TACS  Pitch 

degree 

r###ch04 

05 

TACS  Roll  Rt 

deg/sec 

r###ch05 

06 

TACS  PitchRt 

deg/sec 

r###ch06 

07 

TACS  Yaw  Rt 

deg/sec 

r###ch07 

08 

TACS  CG  ZAcc 

G 

r###ch08 

09 

TACS  CG  YAcc 

G 

r###ch09 

10 

TACS  CG  XAcc 

G 

r###chl0 

11 

TACS  CT  ZAcc 

G 

r###chl  1 

12 

TACS  CT  YAcc 

G 

r###ch!2 

13 

Lhtr  Roll 

degree 

r###chl3 

14 

Lhtr  Pitch 

degree 

r###chl4 

15 

Lntr  Roll  Rt 

deg/sec 

r###chl5 

16 

Lhtr  PitchRt 

deg/sec 

r###ch!6 

17 

Lhtr  CG  ZAcc 

G 

r###chl  7 

18 

Lhtr  CG  YAcc 

G 

r###chl8 

19 

Lhtr  CG  XAcc 

G 

r###ch!9 

20 

TACS  LhtReIX 

inch 

r###ch20 

21 

TACS  LhtRelY 

inch 

r###ch21 

22 

TACS  LhtRelZ 

inch 

r###ch22 

23 

Load  Mo  Orthognl 

inch 

r###ch23 

24 

Load  Mo  Parallel 

inch 

r###ch24 

25 

BoomTip-HorizAcc 

G 

r###ch25 

26 

BoomTip-Vert  Acc 

G 

r###ch26 

27 

Lghtr4-Strn  VAcc 

G 

r###ch27 

28 

Lghtr4-Bow  VAcc 

G 

r###ch28 

29 

Lghtr3-Bow  VAcc 

G 

r###ch29 

30 

Lghtr2-Bow  VAcc 

G 

r###ch30 

31 

Lghtrl-Bow  VAcc 

G 

r###ch31 

32 

Lghtr2-SMdShpVAc 

G 

r###ch32 

33 

Lghtr2-PMdShpVAc 

G 

r###ch33 

34 

Lghtr2-CLMdSpVAc 

G 

r###ch34 

35 

Lghtrl  -CLMdSpVAc 

G 

r###ch35 

36 

Lghtr4-StrnLngAc 

G 

r###ch36 

37 

Lghtr4-PStrnTAcc 

G 

r###ch37 

38 

Lghtrl  -PBowTvAcc 

G 

r###ch38 
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APPENDIX  F.  MASTER  CHANNEL  LIST:  CONFIGURATION  in 


DTMB  Runs  534-587  and  NRL  Runs  227-251 


Channel  # 

Channel  Title 

Units 

Final  File  Title 

01 

Wave  Ht  Bow 

inch 

r###ch01 

02 

Sonix  Sonic 

inch 

r###ch02 

03 

TACS  Roll 

degree 

r###ch03 

04 

TACS  Pitch 

degree 

r###ch04 

05 

TACS  Roll  Rt 

deg/sec 

r###ch05 

06 

TACS  Pitch  Rt 

deg/sec 

r###ch06 

07 

TACS  Yaw  Rt 

deg/sec 

r###ch07 

08 

TACS  CG  ZAcc 

G 

r###ch08 

09 

TACS  CG  YAcc 

G 

r###ch09 

10 

TACS  CG  XAcc 

G 

r###chl0 

11 

TACS  CT  ZAcc 

G 

r###chl  1 

12 

TACS  CT  YAcc 

G 

r###ch!2 

13 

DD  963  Roll 

degree 

r###ch!3 

14 

DD  963  Pitch 

degree 

r###ch!4 

15 

DD  963  CG  ZAcc 

G 

r###chl5 

16 

DD  963  CG  YAcc 

G 

r###ch!6 

17 

DD  963  CG  XAcc 

G 

r###ch!7 

18 

TACS1  DD963  ReIX 

inch 

r###ch!8 

19 

TACS!  DD963  RelY 

inch 

r###chl9 

20 

TACS1  DD963  RelZ 

Inch 

r###ch20 

21 

Load  Mo  Orthognl 

Inch 

r###ch21 

22 

Load  Mo  Parallel 

Inch 

r###ch22 

23 

BoomTip-HorizAcc 

G 

r###ch23 

24 

BoomTip-Vert  Acc 

G 

r###ch24 
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APPENDIX  G.  RAO  PROGRAM 


function  rao(run) 

%  Program  Title:  rao.m - ; - 

%  Referred  to  as  RAO  Program 

%  Must  enter  one  of  the  following  to  initiate  this  program: 
% 

%  rao(##)  for  DTMB  run  numbers  <100 

%  rao(###)  for  DTMB  run  numbers  >100 


%  Initial  Set-up - 

sread  =  char('H:\'); 

swri te  =  char ( ' D : \ ' ) ; 


format  short  e 


ptavg 

window 

lambda 

freq 


50; 

500; 

24.175; 

32.2; 


% 

if 


Establish  Proper  String  Configuration  for  Data  Run 


run  <=  99 

si 

=  char ( ' R0 ' )  ; 

s2 

=  char ( ' 0 ' ) ; 

else 

sl 

=  char ( ' R ' ) ; 

end 


srun 

filename 


num2str  (run)  ; 
s treat (si, srun) ; 


if  run  <=  99 
folder 
else 

folder 

end 


streat (s2 , srun) ; 
srun; 
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%  Determine  Number  of  Data  Channels  to  be  Processed 
if  run  <=503 

channels  =  46; 
elseif  run  <=  533 
channels  =  38; 

else 

channels  =  24; 

end 


Establish  Filepath  and  Load  the  Raw, 


f ilepath 
load (filepath) 
szdata 
channels 


n 


=  [sread, filename] ; 

=  size (data); 

=  szdata (2) ; 

=  szdata (1); 


Time-based  Data 


%  Compute  both  Model-Scale  and  Full-Scale  Omega  Ranges 
momega  =  2*pi*freq/n* (0 :window-l) ; 

fomega  =  momega /sqrt (lambda); 


%  Full-Scale 

f  si 

pul 

rul 


Scalers,  PSD  Units  and  RAO  Units  Library 
=  1; 

=  ' (in) ' ; 

=  ' (dimensionless)  '  ; 


fs2 

pu2 

ru2 


12 /lambda ; 

' (deg)  ' ; 

'  (deg/ in)  '  ; 


f  s3 
pu3 
ru3 


12/ lambda ^ (3/2) ; 
'  (deg/s)  '  ; 

'  (deg/ in  s)  '  ; 


fs4 

pu4 

ru4 


12/lambda^2 ; 

'  (deg/sA2 )  '  ; 

' (deg/in  s^2)  '  ; 


%  Looping/ Processing  of  Data  Channels 
for  ch  =  1: channels 

if  channels==46 


if 


ch==l 

fscaler  = 
psdunits  = 
raounits  = 


fsl; 

pul; 

rul; 
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chtitle  =  'Wave  Ht  Bow'  ; 
elseif  ch==2 

f scaler  =  fsl; 


psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==3 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==4 

fscaler  = 
■psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==5 

fscaler  = 
psdunits  = 
raounits  = 
chtitle 
elseif  ch==6 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==7 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==8 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==9 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==10 
fscaler 
psdunits 
raounits 
chtitle 


pul; 
rul; 

'  Sonix  Sonic ' ; 

fs2; 
pu2 ; 
ru2  ; 

'TACS  Roll'; 

fs2; 
pu2 ; 
ru2 ; 

'TACS  Pitch'; 

fs3; 
pu3 ; 
ru3 ; 

'TACS  Roll  Rt'; 

fs3; 
pu3  ; 
ru3 ; 

'TACS  Pitch  Rt'; 

f  s3  ; 
pu3 ; 
ru3 ; 

' TACS  Yaw  Rt ' ; 

fs4; 
pu4 ; 
ru4; 

' TACS  CG  ZAcc ' ; 

f  s4  ; 

=  pu4 ; 

=  ru4 ; 

=  ' TACS  CG  YAcc ' ; 

=  f  s4 ; 

=  pu4 ; 

=  ru4 ; 

=  ' TACS  CG  XAcc ' ; 
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elseif  ch==ll 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==12 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==13 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==14 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==15 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==16 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==17 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==18 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 
elseif  ch==19 

fscaler  = 
psdunits  = 
raounits  = 
chtitle  = 


fs4; 

pu4; 

ru4; 

'TACS  CT  ZAcc ' ; 

fs4; 
pu4; 
ru4 ; 

' TACS  CF  YAcc ' ; 

fs2; 
pu2; 
ru2  ; 

'Cntr  Roll'; 

fs2; 
pu2 ; 
ru2 ; 

'Cntr  Pitch'; 

fs4; 
pu4  ; 
ru4  ; 

' Cntr  CG  ZAcc ' ; 

fs4; 
pu4 ; 
ru4  ; 

' Cntr  CG  YAcc ' ; 

fs4; 
pu4  ; 
ru4 ; 

' Cntr  CG  XAcc ' ; 

fsl; 

pul; 

rul; 

'TACS  Cntr  Rel  X' 

fsl; 
pul; 
rul ; 

'TACS  Cntr  Rel  Y' 
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elseif  ch==20 

fscaler 

= 

f  si  ; 

psdunits 

= 

pul; 

raounits 

= 

rul  ; 

chtitle 

= 

'TACS 

Cntr  Rel 

elseif  ch==21 

fscaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2  ; 

chtitle 

= 

'  Lhtr 

Roll ' ; 

elseif  ch==22 

fscaler 

= 

f  s2  ; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2 ; 

chtitle 

= 

'Lhtr 

Pitch' ; 

elseif  ch==23 

fscaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'Lhtr 

Roll  Rt' ; 

elseif  ch==24 

fscaler 

= 

fs3; 

psdunits 

- 

pu3  ; 

raounits 

= 

ru3 ; 

chtitle 

= 

'Lhtr 

Pitch  Rt' 

elseif  ch==25 

fscaler 

= 

f  s4 ; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4  ; 

chtitle 

= 

'Lhtr 

CG  ZAcc ' ; 

elseif  ch==2 6 

fscaler 

=s 

fs4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4  ; 

chtitle 

= 

'Lhtr 

CG  YAcc ' ; 

elseif  ch==27 

fscaler 

= 

f  s4 ; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4 ; 

chtitle 

= 

'Lhtr 

CG  XAcc ' ; 

elseif  ch==28 

fscaler 

= 

f  si ; 

psdunits 

pul; 

raounits 

= 

rul; 

chtitle 

= 

'TACS 

Lhtr  Rel 
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elseif  ch==29 

f scaler  =  fsl; 
psdunits  =  pul; 
raounits  =  rul; 
chtitle  =  'TACS  Lhtr  Rel  Y' ; 
elseif  ch==30 

f scaler  =  fsl; 
psdunits  =  pul; 
raounits  =  rul; 
chtitle  =  ' TACS  Lhtr  Rel  Z '  ; 

elseif  ch==31 

f scaler  =  fsl; 
psdunits  =  pul; 
raounits  =  rul; 
chtitle  =  'Load  Mo  Orhognl'; 
elseif  ch==32 

f scaler  =  fsl; 
psdunits  =  pul; 
raounits  =  rul  ; 
chtitle  =  'Load  Mo  Parallel'; 
elseif  ch==33 

fscaler  =  fs4; 
psdunits  =  pu4; 
raounits  =  ru4; 

chtitle  =  'Boom  Tip-HorizAcc ' ; 
elseif  ch==34 

fscaler  =  fs4; 
psdunits  =  pu4  ; 
raounits  =  ru4; 
chtitle  =  'Boom  Tip-VertAcc '  ; 
elseif  ch==35 

fscaler  =  fs4; 
psdunits  =  pu4; 
raounits  =  ru4  ; 
chtitle  =  'Lghtr4-Strn  VAcc ' ; 
elseif  ch==36 

fscaler  =  fs4; 
psdunits  =  pu4; 
raounits  =  ru4 ; 
chtitle  =  'Lghtr4-Bow  VAcc'; 
elseif  ch==37 

fscaler  =  fs4; 
psdunits  =  pu4  ; 
raounits  =  ru4; 
chtitle  =  'Lghtr3-Bow  VAcc'; 
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elseif  ch==38 


f scaler 

= 

fs4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4  ; 

chtitle 

= 

' Lghtr2 -Bow  VAcc ' ; 

elseif  ch==39 

fscaler 

= 

f  s4  ; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4 ; 

chtitle 

= 

'Lghtrl-Bow  VAcc'; 

elseif  ch==40 

fscaler 

= 

fs4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4  ; 

chtitle 

= 

' Lghtr2-SMdShp  VAcc'; 

elseif  ch==41 

fscaler 

= 

fs4  ; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4; 

chtitle 

= 

' Lghtr2-PMdShp  VAcc'; 

elseif  ch==42 

fscaler 

= 

fs4; 

psdunits 

= 

pu4 ; 

raounits 

= 

ru4  ; 

chtitle 

- 

' Lghtr2 -CLMdSp  VAcc ' ; 

elseif  ch==43 

fscaler 

= 

fs4; 

psdunits 

= 

pu4 ; 

raounits 

= 

ru4  ; 

chtitle 

= 

' Lghtrl-CLMdSp  VAcc'; 

elseif  ch==44 

fscaler 

= 

f  s4 ; 

psdunits 

= 

pu4 ; 

raounits 

“ 

ru4 ; 

chtitle 

= 

' Lghtr4-StrnLng  Acc'; 

elseif  ch==45 

fscaler 

= 

fs4; 

psdunits 

= 

pu4 ; 

raounits 

= 

ru4  ; 

chtitle 

= 

' Lghtr4-PStrn  TvAcc'; 

else 

fscaler 

= 

fs4; 

psdunits 

= 

pu4 ; 

raounits 

= 

ru4 ; 

chtitle 

= 

'Lghtrl-PBow  TvAcc'; 

end  %  End  of 

Configuration  I  sub-Loop 

elseif  channels==38 


if  ch==l 


fscaler 

= 

fsl; 

psdunits 

- 

pul; 

raounits 

= 

rul  ; 

chtitle 

= 

'Wave 

Ht  Bow' ; 

elseif  ch==2 

fscaler 

= 

fsl  ; 

psdunits 

= 

pul; 

raounits 

= 

rul; 

chtitle 

= 

' Sonix  Sonic ' ; 

elseif  ch==3 

fscaler 

= 

fs2; 

psdunits 

- 

pu2  ; 

raounits 

= 

ru2  ; 

chtitle 

= 

'TACS 

Roll' ; 

elseif  ch==4 

fscaler 

= 

fs2; 

psdunits 

= 

pu2  ; 

raounits 

= 

ru2 ; 

chtitle 

= 

'TACS 

Pitch' ; 

elseif  ch==5 

fscaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Roll  Rt' 

elseif  ch==6 

fscaler 

= 

fs3; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Pitch  Rt 

elseif  ch==7 

fscaler 

= 

fs3  ; 

psdunits 

= 

pu3  ; 

raounits 

= 

ru3  ; 

chtitle 

= 

'TACS 

Yaw  Rt ' ; 

elseif  ch==8 

fscaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'TACS 

CG  ZAcc ' 
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elseif  ch==9 
f scaler 
psdunits 
raounits 
chtitle 
elseif  ch==10 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==ll 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==12 
fsscaler 
psdunits 
raounits 
chtitle 
elseif  ch==13 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==14 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==15 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==16 
fscaler 
psdunits 
raounits 
chtitle 
elseif  ch==17 
fscaler 
psdunits 
raounits 
chtitle 


fs4; 
pu4 ; 
ru4; 

'TACS  CG  YAcc ' 

f  s4; 
pu4; 
ru4; 

'TACS  CG  XAcc ' 

fs4; 

pu4; 

ru4; 

'TACS  CT  ZAcc ' 

=  f  s4 ; 
pu4  ; 
ru4  ; 

'TACS  CF  YAcc' 

fs2; 
pu2; 
ru2 ; 

'Lhtr  Roll'; 

fs2; 
pu2; 
ru2 ; 

'Lhtr  Pitch' ; 

f  s3 ; 
pu3 ; 
ru3 ; 

'Lhtr  Roll  Rt' 

f  s3  ; 
pu3  ; 
ru3 ; 

'Lhtr  Pitch  RT 

fs4; 
pu4; 
ru4  ; 

'Lhtr  CG  ZAcc' 
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elseif  ch==18 


f scaler 

- 

f  s4 ; 

psdunits 

= 

pu4; 

raounits 

= 

ru4  ; 

chtitle 

= 

'  Lhtr 

CG  YAcc ' ; 

elseif  ch==19 

f scaler 

= 

fs4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4  ; 

chtitle 

= 

'Lhtr 

CG  XAcc ' ; 

elseif  ch==20 

fscaler 

= 

f  si ; 

psdunits 

= 

pul; 

raounits 

= 

rul ; 

chtitle 

= 

'TACS 

LhtRelX' ; 

elseif  ch==21 

fscaler 

= 

fsl; 

psdunits 

= 

pul; 

raounits 

= 

rul; 

chtitle 

= 

'TACS 

LhtRelY' ; 

elseif  ch==22 

fscaler 

= 

fsl ; 

psdunits 

= 

pul; 

raounits 

= 

rul  ; 

chtitle 

= 

'TACS 

LhtRelZ ' ; 

elseif  ch==23 

fscaler 

- 

fsl; 

psdunits 

= 

pul; 

raounits 

= 

rul; 

chtitle 

= 

'  Load 

Mo  Orthognl ' ; 

elseif  ch==24 

fscaler 

= 

fsl; 

psdunits 

pul; 

raounits 

= 

rul ; 

chtitle 

= 

'  Load 

Mo  Parallel'; 

elseif  ch==25 

fscaler 

= 

f  s4  ; 

psdunits 

= 

pu4 ; 

raounits 

= 

ru4 ; 

chtitle 
elseif  ch==26 

— 

' BoomTip-HorizAcc ' ; 

fscaler 

= 

fs4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4 ; 

chtitle 

= 

' BoomTip-Vert  Acc ' ; 
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elseif  ch==27 


fscaler 

= 

f  s4; 

psdunits 

= 

pu4; 

raounits 

= 

ru4; 

chtitle 

= 

'Lghtr4-Strn  VAcc'; 

elseif  ch==28 

fscaler 

= 

fs4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4 ; 

chtitle 

= 

'Lghtr4-Bow  VAcc' ; 

elseif  ch==29 

fscaler 

= 

fs4; 

psdunits 

= 

pu4 ; 

raounits 

= 

ru4; 

chtitle 

= 

' Lghtr3 -Bow  VAcc ' ; 

elseif  ch==30 

fscaler 

= 

fs4  ; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4; 

chtitle 

= 

'Lghtr2-Bow  VAcc'; 

elseif  ch==31 

fscaler 

= 

fs4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4; 

chtitle 

= 

'Lghtrl-Bow  VAcc' ; 

elseif  ch==32 

fscaler 

= 

f  s4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4  ; 

chtitle 

= 

' Lght r 2 - SMdShpVAc ' ; 

elseif  ch==33 

fscaler 

= 

fs4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4  ; 

chtitle 

= 

' Lght r 2 - PMdShpVAc ' ; 

elseif  ch==34 

fscaler 

= 

f  s4 ; 

psdunits 

= 

pu4 ; 

raounits 

= 

ru4  ; 

chtitle 

= 

' Lght r 2 -CLMdSpVAc ' ; 

elseif  ch==35 

fscaler 

= 

fs4; 

psdunits 

= 

pu4  ; 

raounits 

= 

ru4; 

chtitle 

= 

' Lghtrl -CLMdSpVAc ' ; 
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elseif  ch==36 

f scaler  =  fs4; 
psdunits  =  pu4  ; 
raounits  =  ru4 ; 
chtitle  =  ' Lghtr4-StrnLngAc ' ; 

elseif  ch==37 

f scaler  =  fs4; 
psdunits  =  pu4; 
raounits  =  ru4; 
chtitle  =  ' Lghtr4-PStrnTAcc ' ; 

else 

f scaler  =  fs4; 
psdunits  =  pu4; 
raounits  =  ru4; 
chtitle  =  ' Lghtrl-PBowTvAcc ' ; 

end  %  End  of  Configuration  II  sub-Loop 

else 


if  ch==l 


fscaler  = 

f  si ; 

psdunits  = 

pul; 

raounits  = 

rul; 

chtitle  = 

.f  ch==2 

'Wave 

Ht  Bow' ; 

fscaler  = 

f  si  ; 

psdunits  = 

pul; 

raounits  = 

rul ; 

chtitle  = 

elseif  ch==3 

' Sonix  Sonic ' ; 

fscaler  = 

fs2; 

psdunits  = 

pu2  ; 

raounits  = 

ru2  ; 

chtitle  = 

elseif  ch==4 

'TACS 

Roll' ; 

fscaler  = 

fs2; 

psdunits  = 

pu2  ; 

raounits  = 

ru2  ; 

chtitle  = 

elseif  ch==5 

'TACS 

Pitch' ; 

fscaler  = 

fs3; 

psdunits  = 

pu3  ; 

raounits  = 

ru3  ; 

chtitle  = 

'TACS 

Roll  Rt' 
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elseif  ch= 

==6 

f scaler 

=  fs3; 

psdunits 

=  pu3 ; 

raounits 

=  ru3 ; 

chtitle 

=  ' TACS 

Pitch  Rt 

elseif  ch= 

--=1 

fscaler 

=  fs3; 

psdunits 

=  pu3 ; 

raounits 

=  ru3 ; 

chtitle 

=  ' TACS 

Yaw  Rt ' ; 

elseif  ch= 

==8 

fscaler 

=  fs4; 

psdunits 

=  pu4 ; 

raounits 

=  ru4 ; 

chtitle 

=  '  TACS 

CG  ZAcc ' 

elseif  ch= 

:=9 

fscaler 

=  fs4; 

psdunits 

=  pu4 ; 

raounits 

=  ru4 ; 

chtitle 

=  '  TACS 

CG  YAcc ' 

elseif  ch= 

=  =  10 

fscaler 

=  fs4; 

psdunits 

=  pu4 ; 

raounits 

=  ru4  ; 

chtitle 

=  '  TACS 

CG  XAcc ' 

elseif  ch= 

:  =  11 

fscaler 

=  fs4; 

psdunits 

=  pu4 ; 

raounits 

=  ru4 ; 

chtitle 

=  '  TACS 

CT  ZAcc' 

elseif  ch= 

■—12 

fscaler 

=  fs4; 

psdunits 

=  pu4 ; 

raounits 

=  ru4 ; 

chtitle 

=  '  TACS 

CF  YAcc' 

elseif  ch= 

:  =  13 

fscaler 

=  f  s2 ; 

psdunits 

=  pu2 ; 

raounits 

=  ru2 ; 

chtitle 

=  'DD  963  Roll'; 

elseif  ch= 

:  =  14 

fscaler 

=  fs2; 

psdunits 

=  pu2 ; 

raounits 

=  ru2 ; 

chtitle 

=  'DD  963  Pitch' 
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elseif  ch= 

=15 

f scaler 

=  fs4; 

psdunits 

=  pu4 ; 

raounits 

=  ru4  ; 

chtitle 

=  'DD  963 

CG  ZAcc ' ; 

elseif  ch= 

=16 

f scaler 

=  fs4; 

psdunits 

=  pu4 ; 

raounits 

=  ru4  ; 

chtitle 

=  'DD  963 

CG  YAcc ' ; 

elseif  ch= 

=17 

f scaler 

=  fs4; 

psdunits 

=  pu4 ; 

raounits 

=  ru4 ; 

chtitle 

=  'DD  963 

CG  XAcc ' ; 

elseif  ch= 

=18 

fscaler 

=  fsl; 

psdunits 

=  pul  ; 

raounits 

=  rul ; 

chtitle 

=  '  TACS1 

DD963  RelX' 

elseif  ch==19 

fscaler 

=  fsl; 

psdunits 

=  pul  ; 

raounits 

=  rul  ; 

chtitle 

=  '  TACS1 

DD963  RelY' 

elseif  ch= 

:=20 

fscaler 

=  fsl; 

psdunits 

=  pul  ; 

raounits 

=  rul  ; 

chtitle 

=  '  TACS1 

DD963  RelZ ' 

elseif  ch= 

:—21 

fscaler 

=  fsl; 

psdunits 

=  pul ; 

raounits 

=  rul ; 

chtitle 

=  ' Load  Mo  Orthognl ' 

elseif  ch= 

:=2  2 

fscaler 

=  fsl; 

psdunits 

=  pul ; 

raounits 

=  rul ; 

chtitle 

=  'Load  Mo  Parallel' 

elseif  ch= 

-=2  3 

fscaler 

=  fs4; 

psdunits 

=  pu4  ; 

raounits 

=  ru4 ; 

chtitle 

=  ' BoomTip-HorizAcc ' 
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else 

fscaler  =  fs4; 

psdunits  =  pu4  ; 
raounits  =  ru4  ; 

chtitle  =  ' BoomTip-Vert  Acc'; 

end  %  End  of  Configuration  III  sub-Loop 

end  %  End  of  the  Overall  Configuration  Loop 

%  Base-lining  the  Raw,  Time-based  Data - 

x  =  data ( : , ch) ; 

x  =  x-mean(x); 

%  FFT  Computation - 

Y  =  f  f  t  (x,  n)  ; 

%  Model-Scale  PSD  Computation - 

PSD  =  Y. *conj (Y) /n; 

%  Moving-window  Averaging  of  Model-Scale  PSDs  for 

smoothing - 

halfptavg  =  ptavg/ 2; 

sum  =  0  ; 

for  k  =  1: ptavg 

sum  =  sum  +  PSD(k,l); 

if  k  >  halfptavg 

psd (k-halfptavg, ch)  =  sum/k; 
else 
end 

end 

pts  =  window  +  halfptavg; 

ptavgplusone  =  ptavg  +  1; 

3  =  1  ; 

for  k  =  ptavgplusone : pts 

sum  =  sum  +  PSD(ptavg+j , 1)  -  PSD(j,l); 

psd (k-half ptavg , ch)  =  sum/ptavg; 

j  =  j  +  1; 

end 

meansqr  =  mean(x.^2); 

psdarea  =  trap (momega, psd ( : , ch) ) ; 

psdscaler  =  meansqr /psdarea ; 

psd ( : , ch)  =  psd ( : , ch) *psdscaler ; 


%  Computation  of  Model-Scale  RAOs - 

X  =  psd ( : , ch) . /psd (  :  ,  1 )  ; 

mrao(:,ch)  =  sqrt(X); 

%  Computation  of  Full-Scale  RAOs - 

frao(:,ch)  =  mrao (:, ch) *f scaler ; 

%  Plot  Model-Scale  PSDs  and  Save  Figure  in  jpg  Format - 

elf 

plot (momega (1, : ) ,psd ( : , ch) ) 
grid 

sch  =  num2str(ch); 

title ([ 'Model-Scale  PSDs  for  DTMB  Run  'srun,'  /  Channel 
' , sch, ' :  ' chtitle] ) 

xlabel ( ' Omega  ( rad/ s ) ' ) 
ylabel(t'PSD  'psdunits]) 

saveas(gcf , [swrite, folder, '\' , filename, 'C' ,sch, 'p' ] , 'jpg 
');  %  D:\###\R###C##p. jpg 

%  Plot  Model -Scale  RAOs  and  Save  Figure  in  jpg  Format - 

elf 

plot (momega ( 1 , : ) , mrao ( : , ch) ) 
grid 

title ([ 'Model-Scale  RAOs  for  DTMB  Run  'srun,'  /  Channel 
',sch, ':  'chtitle]) 
xlabel ( ' Omega  ( rad/ s ) ' ) 
ylabel(['RAO  'raounits]) 

saveas (gcf , [swrite, folder, ' \ ' , filename, 'C' , sch, 'm' ] , ' jpg 
');  %  D: \###\R###C##m. jpg 

%  Plot  Full-Scale  RAOs  and  Save  Figure  in  jpg  Format - 

elf 

plo t ( f omega ( 1 , : ) , f rao ( : , ch ) ) 
grid 

title ([ 'Full-Scale  RAOs  for  DTMB  Run  'srun,'  /  Channel 
' , sch, ' :  ' chtitle] ) 

xlabel (' Omega  (rad/s)') 
ylabel(['RAO  'raounits]) 

saveas (gcf, [swrite, folder, ' \ ' , filename, ' C ' , sch, ' f ' ] , 'jpg 
');  %  D:\###\R###C##f . jpg 

end  %  End  of  all  Channel  Looping 
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%  Save  Omega  Matrices  as  Last  Column  of  Composite  PSD  and 

RAO  Matrices - - - 

ch  =  ch  +  1; 

psd ( : , ch)  =  momega(l,  : )  '  ; 

mrao ( : , ch)  =  momega ( 1 ,  : )  '  ; 

f rao ( : , ch)  =  f omega (1, : )  '  ; 

%  Save  Composite  PSD  and  RAO  Matrices - 

%  D: \###\R###psd. txt 

%  D: \###\R###mrao . txt 

%  D: \###\R###frao . txt 

save ( [swrite, folder, ' \ ' , filename, 'psd. txt ' ] , 'psd' , ' -ascii 
' , ' -tabs ' )  ; 

save ( [swrite, folder, ' \ ' , filename, 'mrao . txt' ] , 'mrao' , ' -ascii 
' , ' -tabs ' ) ; 

save ( [swrite, folder, ' \ ' , filename, ' frao.txt' ] , ' frao' , ' -ascii 
' , ' -tabs ' ) ; 

%  End  of  RAO  Program - 
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APPENDIX  H.  PTAVG  COMPARISON  PROGRAM 


%  Program  title:  ptavgcomparison . m - 

clear 

clc 

%  Load  Spreadsheet  Results  for  run  67 - 

load  run67 . txt ; 

%  Load  RAO  Program  Results  for  ptavg  =26,  50,  76  and  100-- 

load  R067psd26.txt; 

load  R067psd50.txt; 

load  R067psd76.txt; 

load  R067psdl00.txt; 

load  R067mrao26.txt; 
load  R067mrao50.txt; 
load  R067mrao76.txt; 
load  R067mraol00.txt; 

load  R067frao26.txt; 
load  R067frao50.txt; 
load  R067frao76.txt; 
load  R0 67 f raolO 0 . txt ; 

%  Plot  Model -Scale  PSD  Comparison - 

figure (1) 

plot (R067psd26 ( : , 47) , R067psd26 ( : , 3 ) , . . . 

R067psd50 ( : , 47 ) ,R067psd50 ( : , 3 ) , . . . 

R067psd7  6 ( : ,47) ,R067psd76 ( : ,3) ,  . . . 

R067psdl00 ( : ,47) ,R067psdl00 ( : ,3) , . . . 
run67  (  :  ,  2 )  ,  run67  ( :  ,  4 )  ) 

title ( 'R067CH3  Model-Scale  TACS  Roll  PSDs  "ptavg" 
Comparison' ) 
xlabel ( ' Omega  ( rad/ s ) ' ) 
ylabel ( ' PSD  (deg) ' ) 
grid 
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%  Plot  Model-Scale  RAO  Comparison - 

figure (2 ) 

plot (R067mrao2  6  ( : , 47) , R067mrao26 ( : , 3 )  ,  .  .  . 
R067mrao50 ( : , 47) , R067mrao50 ( : , 3)  ,  .  .  . 
R067mrao76 ( : , 47) , R067mrao76 ( : , 3) , . . . 
R067mraol00 ( : , 47 ) , R067mraol00 (:  ,3)  ,  . .  . 
run67 (: ,2) , run67 ( : , 5 ) ) 

title ( 'R067CH3  Model-Scale  TACS  Roll  RAOs  "ptavg" 
Comparison' ) 
xlabel (' Omega  (rad/s)') 
ylabel('RAO  (deg/in)') 
grid 

%  Plot  Full-Scale  RAO  Comparison - 

figure (3 ) 

plot (R067frao26 ( : , 47) , R067frao26 ( : ,  3 )  ;  .  .  . 
R067frao50 ( : , 47) ;R067frao50 ( : ,3)  ,  .  .  . 
R067frao76 ( : , 47) ,R067frao76 ( : ,3) , . . . 
R067fraol00 ( : ,47) ,R067fraol00 ( :  ,3)  ) 
title ( 'R067CH3  Full-Scale  TACS  Roll  RAOs  "ptavg" 
Comparison' ) 
xlabel (' Omega  (rad/s)') 
ylabel('RAO  (deg/in)') 
grid 

%  End  of  ptavgcomparison  program - 
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APPENDIX  I.  PSD  AND  RAO  COMPARISONS  PROGRAM 


%  Program  Title:  psdandraocomparisons .m - 

%  Computational  and  Graphical  Comparison  between  the 

various  Methods - 

clear 

clc 

format  short 

%  Load  Applicable  Matrices - 

load  R067;  %  Raw,  Time-based  Data 

load  run67.txt;  %  Spreadsheet  PSD  and  RAO  Results 

load  R067psd.txt;  %  RAO  Program  Model-Scale  PSD  Results 

load  R067mrao.txt;  %  RAO  Program  Model-Scale  RAO  Results 

load  R067frao.txt;  %  RAO  Program  Full-Scale  RAO  Results 

%  Base-line  the  Raw  Time-based - - - 

x  =  data ( : , 1 )  ;  %  Wave  Ht  Data 

x  =  x  -  mean(x); 

z  =  data ( : , 3 )  ;  %  TACS  Roll  Data 

z  =  z  -  mean(z)  ; 

%  Computation  of  Wave  Ht  PSDs  and  Associated  Frequencies - 

[hwelch,  fwelch]  =  pwelch(x,  [],[],  8096,32 .2)  ; 

[hperio,  fperio]  =  periodogram (x, [ ] , 8096 , 32 . 2 ) ; 

[hmulti ,  fmulti ]  =  pmtm(x,  [  ]  ,  8096 , 32 . 2)  ; 

meansqrx  =  mean (x. ^2);  %  Wave  Ht  Data 

welchareah  =  trap  (wwelch,  hwelch)  ;  %  Wave  Ht 

perioareah  =  trap (wperio, hperio)  ; 

multiareah  =  trap (wmulti, hmulti) ; 

hwelchscaled  =  hwelch*meansqrx/welchareah; 

hperioscaled  =  hperio*meansqrx/perioareah; 

hmultiscaled  =  hmulti*meansqrx/multiareah; 


%  Computation  of  TACS  Roll  PSDs  and  Associated  Frequencies- 
[rwelch, fwelch]  =  pwelch ( z ,[],[], 8096 , 32 . 2 ) ; 

[rperio, fperio]  =  periodogram ( z , [ ] , 8096 , 32 . 2 ) ; 

[rmulti, fmulti]  =  pmtm(z ,  [  ] , 8096 , 32 . 2 ) ; 

meansqrz  =  mean (z . A2 ) ;  %  TACS  Roll  Data 

welcharear  =  trap (wwelch, rwelch) ;  %  TACS  Roll 

perioarear  =  trap (wperio, rperio) ; 

multiarear  =  trap (wmulti, rmulti) ; 

rwelchscaled  =  rwelch*meansqrz /welcharear; 

rperioscaled  =  rperio*meansqrz /perioarear ; 

rmultiscaled  =  rmulti *meansqrz /multiarear; 

%  Convert  Frequency  to  Omega - 

wwelch  =  2 *pi* fwelch; 

wperio  =  2 *pi* fperio; 

wmulti  =  2 *pi* fmulti; 

%  Compute  Areas  under  the  Wave  Ht  PSD  Curves - 

welchareah2  =  trap (wwelch, hwelchscaled) ; 

perioareah2  =  trap (wperio, hperioscaled)  ; 

multiareah2  =  trap (wmulti, hmultiscaled)  ; 

spreadareah  =  trap (run67 ( : , 2 ) , run67 ( : , 3 ) ) ; 

raoprogareah  =  trap (R067psd ( : ,47) ,R067psd( : , 1) ) ; 

%  Compare  Areas  under  the  Wave  Ht  PSD  Curves  against 

meansqrs - 

wavehtpsdareacomparison  =  [meansqrx  welchareah2  . . . 

perioareah2  multiareah2  spreadareah  raoprogareah] 

%  Compute  Areas  under  the  TACS  Roll  PSD  Curves - - - 

welcharear2  =  trap (wwelch, rwelchscaled)  ; 

perioarear2  =  trap  (wperio,  rperioscaled)  ,- 

multiarear2  =  trap (wmulti,  rmultiscaled) ; 

spreadarear  =  trap (run67 ( : , 2 ) , run67 ( : , 4) ) ; 

raoprogarear  =  trap (R067psd ( : , 47 ) , R067psd(  : , 3 ) ) ; 

%  Compare  Areas  under  the  TACS  Roll  PSD  Curves  against 

meansqrs - 

rollpsdareacomparison  =  [meansqrz  welcharear2  . . . 

perioarear2  multiarear2  spreadarear  raoprogarear] 
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%  Graphical  Model -Scale  Wave  Ht  PSDs  Comparison  of  All  5 
Methods - 

figure (1) 

plot (wwelch, hwelchscaled)  %  Welch  Method 

hold  on 

plot (wperio,hperioscaled)  %  Periodogram  Method 

plot (wmulti,hmultiscaled)  %  Multi-taper  Method 

plot (run67 ( : , 2 ) , run67 ( : ,  3 )  )  %  Spreadsheet  Results 

plot (R067psd ( : , 47 ) , R067psd( : , 1) ) ;  %  Rao  Program  Results 

title ( 'R067CH01  Model-Scale  Wave  Ht  PSDs  Comparison') 

xlabel ( 'Omega  (rad/s) ' ) 

ylabel ( 'PSD  (dimesionless ) ' ) 

axis ([060  . 015] ) 

grid 

%  Graphical  Model-Scale  Wave  Ht  PSDs  Comparison - 

%  Spreadsheet  vs.  RAO  Program  only- 
figure  (2) 

plot (run67 ( : , 2) , run67 ( : , 3) )  %  Spreadsheet  Results 

hold  on 

plot (R067psd( : , 47) ,R067psd( : , 1) )  %  Rao  Program  Results 

title ( 'R067CH01  Model-Scale  Wave  Ht  PSDs  Comparison') 

xlabel (' Omega  (rad/s)') 

ylabel ( 'PSD  (dimesionless) ' ) 

axis ([060  . 0065] ) 

grid 

%  Graphical  Model-Scale  TACS  Roll  PSDs  Comparison  of  All  5 

Methods - 

figure (3 ) 

plot (wwelch, rwelchscaled)  %  Welch  Method 

hold  on 

plot (wperio, rperioscaled)  %  Periodogram  Method 

plot (wmulti,rmultiscaled)  %  Multi-taper  Method 

plot (run67 ( : , 2) , run67 ( : , 4) )  %  Spreadsheet  Results 

plot (R067psd( : ,47) ,R067psd( : ,3) ) ;  %  Rao  Program  Results 

title ( 'R067CH03  Model-Scale  TACS  Roll  PSDs  Comparison') 

xlabel ( 'Omega  (rad/s) ' ) 

ylabel ('PSD  (deg)') 

axis ( [0  6  0  . 02] ) 

grid 
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%  Graphical  TACS  Roll  PSDs  Comparison - 

%  Spreadsheet  vs .  RAO  Program  only 
figure (4) 

plot (run67 ( : , 2) , run67 ( : , 4) )  %  Spreadsheet  Results 

hold  on 

plot (R067psd( : ,47) , R067psd ( : , 3 ) )  %  Rao  Program  Results 

title (' R067CH03  Model-Scale  TACS  Roll  PSDs  Comparison') 

xlabel (' Omega  (rad/s)') 

ylabel('PSD  (deg)') 

axis ([060  . 009] ) 

grid 


%  Compute 
Rl 

welchmrao 

R2 

periomrao 

R3 

multimrao 


Model-Scale  Roll  RAOs  from  MATLAB  Methods 
=  rwelchscaled. /hwelchscaled; 

=  sqrt (Rl) ; 

=  rperioscaled. /hperioscaled; 

=  sqrt (R2 ) ; 

=  rmultiscaled. /hmultiscaled; 

=  sqrt(R3); 


%  Graphical  Model-Scale  TACS  Roll  RAOs  Comparison  of  All  5 

Methods - . - 

figure ( 5 ) 

plot (wwelch, welchmrao) 
hold  on 

plot (wperio , periomrao ) 
plot (wmulti, multimrao) 
plot (run67 ( : , 2 ) , run67 ( : , 5 ) ) 


%  Welch  Method 

%  Periodogram  Method 
%  Multi-taper  Method 
Spreadsheet  Results 


plot (R067mrao ( : , 47 ) , R067mrao ( : , 3 ) )  ;  %  Rao  Program  Results 


title (' R067CH03  Model-Scale  TACS  Roll  RAOs  Comparison') 

xlabel ( 'Omega  (rad/s) ' ) 

ylabel('RAO  (deg/in)') 

axis ( [0  6  0  10] ) 

grid 

%  Graphical  Model-Scale  TACS  Roll  RAOs  Comparison - 

%  (Spreadsheet  vs.  RAO  Program  only) 
figure ( 6 ) 

plot (run67 ( : , 2 ) , run67 ( : , 5 ) )  %  Spreadsheet  Results 

hold  on 

plot (R067mrao ( : , 47 ) , R067mrao ( :  ,  3  )  )  %  Rao  Program  Results 
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title ( 'R067CH03  Model-Scale  TACS  Roll  RAOs  Comparison') 

xlabel ( ' Omega  ( rad/ s ) ' ) 

ylabel ( 'RAO  (deg/ in)  '  ) 

axis ( [ 0  6  0  2.2]) 

grid 


%  Scale  to 
welchfrao 
periofrao 
multifrao 


Full-Scale  RAOs  and  Omegas - 

=  welchmrao*12/24 . 175; 
=  periomrao*12/24 . 175 ; 
=  mul timrao* 12/24.175; 


wwelchfull 

wperiofull 

wmultifull 


wwelch/sqrt (24 . 175)  ; 
wperio/sqrt (24 . 175 )  ; 
wmulti/sqrt (24 . 175)  ; 


%  Graphical  Full-Scale  TACS  Roll  RAOs  Comparison - - 

%  (4  Methods  Not  including  Spreadsheet) 

figure (7) 

plot (wwelchfull, welchfrao)  %  Welch  Method 


hold  on 

plot (wperiofull 
plot (wmultifull 
plot (R067frao ( : 


periofrao) 

multifrao) 

47) ,R067frao(: ,3) ) 


%  Periodogram  Method 
%  Multi- taper  Method 
%  RAO  Program  Results 


title ( 'R067CH03  Full-Scale  TACS  Roll  RAOs  Comparison') 

xlabel ( 'Omega  (rad/s) ' ) 

ylabel ( ' RAO  (deg/in)') 

axis ( [0  1.4  0  3] ) 

grid 


%  Graphical  Full-Scale  TACS  Roll  RAOs  Presentation - 

%  RAO  Program  only 
figure (8) 

plot (R067frao( : ,47) ,R067frao( :  ,3)  )  %  Rao  Program  Results 

title ( 'R067CH03  Full-Scale  TACS  Roll  RAOs  Presentation') 

xlabel ( 'Omega  (rad/s) ' ) 

ylabel ( 'RAO  (deg/in) ' ) 

axis ([0  1.4  0  1.4]) 

grid 


%  End  of  PSD  and  RAO  Comparison  Program' 
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